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Executive  Summary  of  "Vibration  Dynamics  and  Control  of  Bladed  Disk  Assemblies” 

This  final  report  documents  the  work  performed  at  Purdue  University  during  the  period  of 
November  1988  to  December  1990.  The  original  AFOSR  Contracts  (#AFOSR- 89-0002, 
AFOSR- 89-00 14)  were  written  for  two  years.  Consequently  this  research  was  partially  funded 
from  Professor  Nwokahs’  PRF  grant  #670-1667.  The  objective  of  the  proposed  research  was  to 
gain  a  fundamental  understanding  of  how  and  why  periodically  configured  mechanical  and 
structural  systems,  (in  particular  bladed-disk  assemblies)  with  cyclic  synunetry  and  nominally 
identical  sub-structures  can  display  non-uniform  amplitudes  of  vibration  when  subjected  to  small 
but  random  parameter  perturbations  that  are  often  within  the  component  manufacturing 
tolerances.  A  secondary  aim  of  the  proposal  was  to  determine  ways  of  passively/actively  (if 
possible)  controlling  these  uneven  vibration  amplitudes.  This  work  specifically  dealt  with  the 
influence  of  the  double  (degenerate)  eigenvalues  present  in  every  cyclic  mechanical  system  and 
their  subsequent  splitting  under  small  perturbations,  on  the  uneven  vibration  amplitudes  of  the 
components. 

Status: 

The  work  associated  with  the  principal  objectives  of  the  project  is  almost  completed  and  is 
included  in  this  final  report.  The  procedure  for  detecting  a  priori  which  degenerate  eigenvalue 
pairs  will  split  under  given  parameter  perturbations  has  been  formalized  by  use  of  finite  group 
representation  theory  and  is  presented  in  Appendix  1.  The  procedure  for  accurately  unfolding 
the  singularities  induced  by  the  splitting  of  the  double  modes  of  cyclic  systems  has  been 
formalized  by  the  use  of  a  singular  perturbation  analysis  technique  which  is  valid  for  any  finite 
order  cyclic  system  and  is  included  in  Appendix  2. 
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The  topological  basis  for  the  singularities  induced  by  the  double  modes  and  the 
consequences  there  of  are  carefully  examined  and  detailed  in  Appendix  3. 

In  contra-distinction  from  recent  work  in  the  bladed-disk  research  literature,  numerical 
studies  which  show  that  uneven  amplitudes  of  vibration  in  perturbed  cyclic  systems  can  arise 
both  under  strong  coupling  as  well  as  the  weak  coupling  conditions  is  included  in  Appendix  4. 

A  systematic  framework  has  now  been  established  for  a  detailed  study  of  perturbed  cyclic 
systems.  Future  efforts  will  be  aimed  at  completing  any  remaining  theoretical  analysis, 
development  of  computational  algorithms  for  such  analysis,  passive  structural  redesign  to  avoid 
localized  high  vibration  amplitudes  and  experimental  validation  of  the  analytical  results. 

Publications 

Six  papers  have  been  developed  from  this  work.  The  first  is  essentially  Appendix  1.  The 
second  and  third  are  included  in  Appendix  2.  The  fourth  is  given  in  Appendix  3,  while  the  rest 
are  in  conference  proceedings  as  given  below. 

1.  Happawana,  G.S.,  Bajaj,  A.K.,  Nwokah,  O.D.I.,  On  the  dynamics  of  permrbed  symmetric 
systems.  Accepted  for  presentation  and  publication  in  conference  proceedings  for  13th 
Biennial  ASME  Conference  on  Mechanical  Vibration  and  Noise,  September  22-25,  1991, 
Miami,  Florida. 

2.  Happawana,  G.S.,  Bajaj,  A.K.,  Nwokah,  O.D.I.,  A  singular  perturbation  perspective  on 
mode  localization.  J.  Sound  and  Vibration,  (To  appear). 

3.  Happawana,  G.S.,  Bajaj,  A.K.,  Nwokah,  O.D.I.,  A  singular  perturbation  analysis  of 
eigenvalue  veering  and  mode  localization  in  linear  chain  and  cyclic  systems. 
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J.  Sound  and  Vibration,  (Submined). 

4.  Nwokah,  O.D.I.,  Afolabi,  D.,  Damra,  F.M.,  On  the  modal  stability  of  imperfect  cyclic 
systems.  Control  and  Dynamic  Systems,  35, 137-164, 1990. 

5.  Afolabi,  D.,  Nwokah,  O.D.I.,  The  frequency  response  of  mistuned  periodic  systems.  Proc. 
12th  ASME  Biennial  Vibrations  Conference,  Montreal,  Canada,  September  1989. 

6.  Afolabi,  D.,  Nwokah,  O.D.I.,  Effects  of  mild  perturbations  on  the  dynamics  of  structures 
with  circulant  matrices.  Proc.  AIAA/ASME/ASCE/AHS  Structures,  Structural  Dynamics 
and  Materials  Conference,  Long  Beach,  CA,  1990. 

Personnel 

Three  faculty  members  and  one  graduate  student  were  funded  by  this  contract. 

Professor  Osita  D.I.  Nwokah  (Resume  at  the  back). 
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Gemunu  S.  Happawana  (Doctoral  Student). 

A  doctorate  degree  is  expected  to  be  awarded  for  this  work  in  the  next  24  months. 

Presentations 

Several  seminar  presentations  resulted  from  this  work.  Details  are  contained  in  the 


individual  Professors’  resumes. 
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1.  INTRODUCTION 

1.1  Problem  Statement 

The  central  aim  of  the  present  project  has  been  to: 

(i)  Gain  a  fundamental  understanding  of  how  and  why  periodically  configured  mechanical 
and  structural  systems  with  cyclic  symmetry  and  nominally  identical  sub-structures  can 
display  non-uniform  amplitudes  of  vibration  under  differential  (i.e.,  small)  parameter 
perturbations  that  are  often  within  the  component  manufacturing  tolerances. 

(ii)  Design  passive  and/or  active  control  mechanisms  to  overcome  such  possible  uneven 
amplitudes  of  vibrations. 

1.2  Background  and  Overview 

The  study  of  cyclically  configured  dynamical  systems,  otherwise  known  as  bladed-disk 
assemblies,  has  been  a  very  active  area  of  research  in  structural  dynamics  over  the  last  25  years. 
It  is  a  measure  of  the  theoretical  difficulties  involved  in  an  accurate  analysis  that  at  the  present 
time  there  is  no  general  agreement  in  the  literature  as  to  either  the  causes  of  the  uneven 
component  vibration  amplitudes  or  as  to  which  component  will  vibrate  with  the  highest 
amplitude  under  parameter  variations.  This  is  a  sine'  qua-non  to  establishing  benchmark 
specifications  for  component  vibration  control.  The  work  performed  under  this  grant  in  the  last 
two  years  clearly  indicates  that: 

(i)  Uneven  amplitudes  of  vibration  are  caused  by  the  modal  bifurcation  phenomenon  or  the 
sensitive  dependence  of  eigenvectors  on  small  parameter  variations  under  some  clearly 
defined  conditions.  [1,2,3] 


-6- 


(ii)  Modal  bifurcations  in  turn  are  caused  by  mistuning  or  small  parameter  variations  from 
nominal  design  values,  which  are  often  within  the  component  manufacturing  tolerances. 
[4] 

(iii)  Extreme  cases  of  the  uneven  amplitudes  of  component  vibrations  produce  the  mode 
localization  phenomenon. 

We  had  shown  in  the  paper  in  Appendix  3,  that  very  useful  qualitative  information  on 
the  blade  mistuning  could  be  obtained  by  application  of  the  methodology  of  singularity 
theory  to  this  problem. 

To  understand  mode  localization,  one  must  first  study  modal  bifurcations. 

Let  |i  e  r  c  R^  where  is  an  r-dimensional  parameter  space.  If  a  given  structural  system 
has  n  degrees  of  freedom,  then  the  characteristic  equation  for  natural  frequencies,  co^  =  can  be 
written  as  the  n-th  order  polynomial  equation: 

G(>, ,  |J.)  =  0 ,  for  some  p.  e  F . 

It  turns  out  that  the  characteristic  polynomial  under  appropriate  modifications  behaves  like  the 
potential  function  in  singularity  (or  catastrophe)  theory.  [5]  Hence  the  degenerate  critical  points 
of  G(X  ,  p)  =  0  correspond  to  the  repeated  roots  (i.e.,  the  repeated  eigenvalues)  of  G(X  ,  p).  By 

studying  together  with  appropriate  higher  order  differentials,  and  G(k ,  p)  =  0,  we  can 
oA, 

determine  the  set  of  all  p  e  F  at  which  G(X ,  p)  has  degenerate  eigenvalues.  This  set,  which  is 
called  the  bifurcation  (or  catastrophe)  set,  partitions  the  parameter  space  F  into  distinct 
submanifolds  whose  boundary  is  the  bifurcation  set. 


We  can  conclude  from  the  basic  theorems  and  results  from  singularity  theory  and  catastrophe 
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theory  [6,7,8]  that  the  modal  behavior  of  our  structural  system  displays  sensitive  dependence  on 
parameters  only  in  the  neighborhood  of  the  bifurcation  set.  We  have  identified  two  distinct 
degenerate  behavior  patterns  in  structural  systems,  namely: 

a)  Coupling  induced  degeneracy, 

b)  Geometry  or  symmetry  induced  degeneracy. 

Furthermore,  we  have  noted  that  one  dimensional  lattice  type  periodic  structures  need  to  be 
divided  into  two  main  classes: 

(i)  The  Linear  Chain, 

(ii)  The  Cyclic  Chain. 

Each  of  these  classes  has  its  own  peculiar  characteristics  which  are  dictated  both  by  the 
geometry  (boundary  conditions)  and  the  physics  of  the  system.  For  example,  in  the  linear  chain, 
degenerate  and  therefore  ‘seemingly’  unpredictable  behavior  under  perturbations  appears  to 
occur  only  under  very  weak  coupling  conditions.  Topologically,  this  behavior  is  equivalent  to 
an  unfolding  of  the  m-fold  (here  m  is  the  number  of  nominally  identical  subsystems  which  are 
weakly  coupled)  degeneracy:  (k  -  O)^)"*  =  0.  This  corresponds  to  what  Pierre  has,  in  a  series  of 
papers,  consistently  referred  to  as  a  perturbation  of  the  uncoupled  system  behavior.  [9,10,11] 
This  behavior  does  not  exist  under  strong  coupling  conditions. 

On  the  other  hand  for  cyclic  systems,  even  under  very  strong  coupling  conditions,  extra 
degeneracy  is  induced  by  the  cyclic  symmetric  nature  of  the  system  matrices.  It  is  then  well 
known  that  cyclic  systems  have  several  pairs  of  degenerate  (coincident)  eigenvalues  which  is 
distinct  from  the  case  of  linear  chains  where  no  degeneracy  or  multiplicity  of  eigenvalues  arises 
under  strong  coupling  conditions.  The  crucial  observation  is  that  because  of  the  coincidence  of 
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eigenvalues,  and  the  continuity  of  eigenvalues  with  respect  to  parameters,  cyclic  systems  (to 
which  bladed-disk  assemblies  belong)  always  operate  in  the  neighborhood  of  the  bifurcation  set. 
For  cyclic  systems  it  is  therefore  of  great  interest  to  determine  the  relative  influence  of  coupling 
and  geometry  in  the  subsequent  degenerate  system  behavior.  Since  a  tuned  bladed-disk  assembly 
has  pairs  of  degenerate  eigenvalues,  the  parameters  corresponding  to  the  tuned  state  are  clearly  a 
subset  of  the  bifurcation  set. 

The  number  of  degenerate  pairs  of  eigenvalues  as  well  as  the  effect  of  different  types  of 
perturbations  depends  on  the  nature  of  the  symmetry.  Some  of  the  qualitative  ramifications  of 
the  geometric  symmetry  can  be  studied  using  the  theory  of  groups.  The  effects  of  symmetry 
preserving  and  symmetry  breaking  perturbations  can  be  qualitatively  studied  using  the  ideas 
from  perturbation  of  group  action  as  well  as  the  singularity  theory  for  symmetric  systems.  While 
the  results  for  universal  unfolding  of  positive  definite  matrices  and  the  behavior  of  eigenvalues 
for  symmetry  preserving  penurbations  are  available  [15],  those  for  penurbations  that  destroy 
symmetry  are  not,  and  we  will  later  present  some  examples  displaying  the  interesting 
consequences  of  various  types  of  perturbations.  Finally,  neither  the  group  theory,  nor  the 
singularity  theory,  provide  quantitative  results  such  as  formulas  for  the  computation  of  the 
perturbed  eigenvalues  and  eigenvectors  as  a  function  of  the  perturbation  parameters.  Only  such 
information  can  provide  the  measures  for  eigenvalue  loci  veering  and  mode  localization,  and  one 
possible  tool  for  developing  these  expressions/results  is  the  singular  perturbation  theory.  Thus, 
tools  or  ideas  from  the  disciplines  of  group  theory,  singularity  theory,  and  singular  perturbation 
theory,  are  all  needed  to  make  a  strong  headway  in  understanding  the  phenomenon  of  mode 


localization. 
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2.  GROUP  THEORY  AND  CYCLIC  SYMMETRY 

Although  it  had  been  observed  that  turned  bladed-disk  assemblies  always  have  many  pairs  of 
degenerate  eigenvalues,  no  theoretical  justification  for  this  phenomenon  was  available  in  the 
bladed-disk  literature.  Our  first  order  of  business  in  the  investigation  was  therefore  to  obtain  a 
formal  explanation  for  this  phenomenon.  The  coefficient  matrices  in  the  equations  of  motion  of 
forced  bladed-disk  assemblies  as  well  as  the  dynamic  stiffness  matrices  are  always  banded 
circulant  matrices  [12].  These  matrices  have  unique  symmetry  properties  [13]  which 
immediately  indicate  that  group  theory  would  be  applicable.  It  turns  out  that  the  set  of  allowed 
symmetry  operations  in  a  bladed-disk  assembly  namely:  rotations  about  a  fixed  axis,  reflections 
about  a  fixed  axis  and  vibrations  about  a  reference  point,  can  be  captured  by  the  operations  of  the 
Dihedral  group  D„  [14].  By  purely  formal  arguments  from  group  theory  and  standard  results  for 
the  Dihedral  group,  we  are  able  to  show  the  number  and  order  of  degenerate  eigenvalues  which 
any  finite  order  bladed-disk  assembly  can  have.  Furthermore  by  considering  the  irreducible 
representations  to  which  the  translational,  rotational,  and  vibrational  modes  belong,  along  with 
the  corresponding  Hamiltonians,  we  can  sufficiently  study  the  effects  of  mild  perturbations  on 
these  degenerate  doublets.  For  example  under  a  given  parameter  variation,  the  symmetry 
operations  generate  a  new  group  which  is  necessarily  a  subgroup  of  the  original  group  Dn.  By 
comparing  the  properties  of  this  new  sub-group  with  those  of  the  original  group,  we  are  able  to 
determine  if  such  a  perturbation  would  lead  to  a  splitting  of  any  of  the  degenerate  pairs  of 
eigenvalues.  It  may  then  be  possible  to  determine  the  minimum  number  of  parameters  which 
must  be  varied  simultaneously  in  order  for  a  certain  number  of  degenerate  eigenvalue  pairs  to  be 
split  at  the  same  time.  By  now  concentrating  on  those  perturbations  that  lead  to  splitting  of 
degenerate  pairs  we  can  more  fully  study  the  effects  of  these  penurbations  on  the  forced 
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amplitude  response  of  the  assembly.  These  results  arc  summarized  in  the  paper  in  Appendix  1. 

3.  SINGULARITY  THEORY  AND  CYCLIC  SYMMETRY 

It  is  known  from  singularity  theory  that  the  splitting  of  eigenvalues  of  a  matrix  can  lead  to 
rapid  changes  in  the  eigenvectors,  which  in  turn  can  result  in  significant  changes  in  the  forced 
amplitude  of  response  of  the  assembly  to  external  aerodynamic  loading.  Since  the  group  theory 
results  indicate  that  bladed-disk  assemblies  could  only  have  degenerate  pairs  of  eigenvalues,  the 
simplest  essential  properties  of  any  finite  order  bladed-disk  assembly  are  captured  by  the 
properties  of  an  assembly  of  order  3.  Note  that  we  need  a  3rd  order  assembly  in  order  to  inscribe 
a  circle  and  hence  obtain  a  cyclic  system.  A  third  order  tuned  assembly  would  thus  have  a 
degenerate  pair  of  eigenvalues  and  an  isolated  eigenvalue.  We  may  therefore  study  the  influence 
of  penurbations  in  the  masses,  ground  springs  and  coupling  springs  on  the  dynamics  of  this 
system.  From  Arnold’s  results  in  singularity  theory  [16],  it  is  self  evident  that  under  mild 
parameter  variations  interest  should  be  concentrated  not  on  the  isolated  eigenvalue  but  only  on 
the  subsequent  behavior  of  the  degenerate  doublet.  To  understand  its  behavior  it  is  necessary  to 
study  the  behavior  of  any  arbitrary  doublet  and  the  subsequent  eigenloci  as  a  function  of 
parameters  in  a  manner  reminiscent  of  root  loci  behavior  in  classical  control  theory.  The 
simplest  doublet  which  contains  the  essential  ingredients  of  the  problem  turns  out  to  be  the 
symmetric,  coupled  double  pendulum  shown  in  figure  1  in  Appendix  2.  The  essense  of  this 
study  was  to  discover  the  relation  of  the  eigenloci  to  parameter  variations  and  the  corresponding 
eigenvectors.  We  had  conjectured  that: 

(i)  Uneven  amplitudes  of  vibration  in  symmetric  structural  systems  are  caused  by  the 
sensitive  dependence  of  system  eigenvectors  or  parameters. 
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(ii)  Sensitive  dependence  of  eigenvectors  on  parameters  gives  rise  under  appropriate 
conditions  to  the  mode  localization  phenomenon. 

(iii)  Rapid  convergence-divergence  (veering)  of  eigenvalues  is  a  signature  for  the  sensitive 
dependence  of  eigenvectors  on  parameters  and  hence  of  the  possible  existence  of  mode 
localization  under  appropriate  conditions. 

If  the  conjecture  were  to  be  true,  we  hoped  to  oe  able  to  obtain  an  estimate  for  the 
eigenvector  sensitivity  measure  in  an  appropriate  manner  as  well  as  an  estimate  for  the 
eigenvector  rotations  resulting  from  any  mild  perturbation.  If  extreme  imperfection  sensitivity 
were  present,  it  was  expected  that  both  measures  would  show  a  singularity  which  is  an  indication 
of  imperfection  sensitivity.  Furthermore  these  were  expected  to  occur  at  the  parameter  values  of 
maximum  curvature  of  the  eigenloci.  If  the  double  pendulum  were  decoupled  (no  coupling 
spring)  there  then  would  exist  two  independent  but  equal  vibration  frequencies.  By  including 
very  weak  coupling  between  the  masses  we  could  study  system  behavior  in  the  neighborhood  of 
the  erstwhile  equal  eigenvalues.  The  study  of  weakly  coupled  systems  is  very  important  since  in 
practice  the  aim  has  always  been  towards  use  of  rigid  disks,  in  effect  making  the  inter-blade 
coupling  very  weak  indeed.  By  assuming  that  the  imperfection  parameter  is  a  slight  difference 
in  the  length  of  the  two  pendula,  we  could  ihei*  ^jdy  the  behavior  of  the  eigenvalues  and 
eigenvectors  of  this  simple  symmetric  system  under  slight  changes  in  coupling  and  disorder. 

4.  QUANTITATIVE  UNFOLDING  OF  THE  MODAL  SINGULARITIES  BY  SINGULAR 
PERTURBATION  ANALYSIS 

We  may  therefore  write  down  the  characteristic  polynomial  as  a  function  of  both  the 
coupling  parameter  and  the  imperfection  parameter.  The  characteristic  polynomial  in  turn 
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behaves  identically  to  a  potential  function  in  singularity  theory  [17].  Thus  the  degenerate 
critical  points  of  this  function  correspond  to  the  repeated  eigenvalues  if  any.  By  writing  down 
the  expressions  for  the  eigenvalues  as  functions  of  the  two  parameters,  the  detailed  behavior  of 
the  eigenloci  in  any  neighborhood  can  be  obtained.  A  regular  expansion  of  these  eigenvalues  as 
a  function  of  the  two  parameters  breaks  down  (loses  uniformity)  in  the  neighborhood  of  the 
critical  point  (where  the  eigenvalues  are  coincident).  By  applying  the  techniques  of  singular 
perturbation  analysis  and  appropriate  stretching  transformations  it  became  possible  to  obtain  the 
eigenloci  expressions  which  were  uniformly  valid  over  the  dv  lain  of  definition  of  the  small 
parameters  and  whose  loci  clearly  indicated  the  veering  phenomenon.  The  same  technique  was 
also  applied  to  the  eigenvector  expressions  as  functions  of  the  parameters.  From  these 
expansions,  expressions  were  obtained  both  for  the  modal  sensitivity  measure  and  the 
eigenvector  rotation  measure  under  slight  parameter  variations.  All  the  results  obtained, 
confirmed  the  conjecture.  The  first  part  of  these  results  are  to  appear  in  the  Journal  of  Sound  and 
Vibration  while  the  second  part  involving  the  full  eigenvector  work  has  been  submitted  to  the 
Journal  of  Sound  and  Vibration.  These  manuscript  preprints  are  enclosed  in  Appendix  2.  We 
can  now  claim  that  we  understand  fairly  well  the  causes  of  localization  phenomenon  for  simple 
doublets.  We  note  however  that  the  double  pendulum  analysis  displayed  the  noted  strange 
behavior  only  under  very  weak  coupling  conditions.  Our  work  had  shown  earlier  (see  Appendix 
4)  that  for  bladed-disk  assemblies  uneven  amplitudes  of  vibration  and  hence  mode  localization 
could  occur  even  under  very  strong  coupling  conditions.  We  therefore  had  to  discover  under 
what  conditions  the  double  pendulum  analysis  remained  valid  also  for  the  bladed-disk  assembly. 
It  turns  out  that  the  pendulum  analysis  remains  valid  for  the  bladed-disk,  irrespective  of  the 
coupling  strength.  However  by  systematically  reducing  the  coupling,  more  complicated 
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singularities  appear.  This  is  because  at  very  low  coupling  the  bladed-disk  behaves  like  a 
penurbation  of  a  triple  degeneracy  which  in  general  requires  at  least  three  parameters  to  unfold 
(completely  analyse).  On  the  other  hand  for  the  double  pendulum,  two  parameters  were  enough 
to  unfold  the  doublet  degeneracy. 

The  lessons  learned  so  far  are  thus  that  for  the  bladed-disks  even  under  strong  coupling 
conditions  mode  localization  can  occur.  On  the  other  hand  for  the  double  pendulum  or  linear 
chains  in  general,  mode  localization  occurs  only  under  weak  coupling.  This  finding  contradicts 
the  current  view  in  bladed-disk  research  [11],  which  bolds  that  in  both  linear  chains  (coupled 
pendula)  and  cyclic  chains  (bladed-disk  assemblies)  mode  localization  only  occurs  under  very 
weak  coupling.  What  is  however  true  for  bladed-disk  assemblies,  is  that  under  very  weak 
coupling  new  singularity  types  (which  do  not  exist  under  strong  coupling)  appear.  We  do  not  yet 
understand  the  full  effects  of  these  new  singularity  types.  We  however  conjecture  that  they  will 
funher  complicate  the  modal  behavior  of  the  assembly  under  aerodynamic  loading.  The  key 
question  we  seek  to  answer  presently  is  which  of  either  symmetry  breaking  perturbations  or 
coupling  induced  perturbations  have  more  influence  on  the  modal  behavior  of  a  bladed-disk 
assembly.  Are  there  regimes  where  each  has  more  influence  than  the  other  and  if  so,  what  is  the 
transition  region?  If  we  could  answer  these  questions  then  we  could  specify  apriori  the 
acceptable  range  of  coupling  so  that  design  effort  could  be  concentrated  on  symmetry  breaking 
bifurcations  and  how  to  prevent  their  effects  from  being  felt  at  the  blade  amplitudes.  The 
singular  perturbation  analysis  acts  as  an  unfolding  of  the  singularities  involved,  since  by  this 
methodology  we  are  able  to  obtain  detailed  information  on  the  modal  behavior  of  the  structure  in 
the  neighborhood  of  the  singularities. 
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5.  THE  CONTROL  PROBLEM 

We  have  not  carried  out  the  control  design  component  of  the  project  as  stated  in  the 
statement  of  work  because  it  has  only  been  in  the  last  six  to  eight  months  that  a  thorough 
understanding  of  the  structural  dynamics  has  emerged.  We  however  are  clear  on  the  work  that 
needs  to  be  done.  The  next  stage  of  our  work  will  involve  classification  of  different 
perturbations  with  the  corresponding  amplitudes  of  vibration.  The  control  problem  in  one 
possible  approach  is  a  structural  redesign  that  deliberately  breaks  the  symmetry  by  splitting  the 
degenerate  eigenvalues  with  only  those  perturbations  that  do  not  lead  to  amplification  of 
vibration  amplitudes.  Provided  the  split  eigenvalues  are  not  in  the  neighborhood  of  the 
bifurcation  set,  all  further  slight  perturbations  would  not  be  expected  to  display  extreme 
imperfection  sensitivity.  Another  alternative  control  methodology  which  we  are  presently 
considering  is  a  regular  adapative  control  scheme  that  seeks  by  means  of  active  addition  or 
subtraction  of  control  masses  and  springs  to  restore  symmetry  whenever  the  symmetry  breaking 
signature  is  observed.  Under  this  scheme  the  degree  of  sensitivity  and  eigenvector  rotation  will 
determine  the  amount  of  modification  called  for  and  the  location  where  to  apply  it.  However 
this  kind  of  scheme  seems  to  us  to  be  more  appropriate  for  aerospace  structural  systems  than  to 
turbine  rotor  disks. 

6.  CONCLUSIONS 

The  results  obtained  from  the  study  in  the  last  two  years  have  helped  to  clarify  and  unify 
several  conflicting  viewpoints  within  the  bladed-disk  research  community.  What  is  more 
significant  is  that  it  has  led  to  a  better  understanding  of  the  potentially  very  complicated 
dynamical  structure  which  ensues  when  geometric  (spatial)  symmetry  interacts  with  weak 
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coupling  in  periodic  structures.  Without  this  understanding  any  attempts  at  either  structural 
redesign  or  structural  control  of  such  systems  would  inevitably  be  frought  with  danger.  We  arc 
currently  continuing  work  on  the  forced  response  of  cyclic  systems  with  a  view  to  a  more 
complete  mathematical  characterization  of  the  relationship  between  amplitudes  of  vibration, 
mode  localization,  and  perturbation  type.  We  are  also  generalizing  the  singular  perturbation 
approach  to  linear  chains  and  cyclic  symmetric  systems  of  any  finite  order.  We  believe  that  the 
development  of  the  structural  control  schemes  would  be  worthless  without  this  full 
understanding. 

REFERENCES 

1.  Afolabi,  D.,  Nwokah,  O.D.I.,  Effects  of  Mild  Perturbations  on  the  Dynamics  of  Structures 
with  Circulant  Matrices,  Proc.  AIAA/ASME/ASCE/AHS  Structures,  Structural  Dynamics 
and  Materials  Conference,  Long  Beach,  CA  1990. 

2.  Afolabi,  D.,  Vibration  Amplitudes  of  Mistuned  Blades,  J.  Turbomachinery,  110,  251-257, 
1988. 

3.  Nwokah,  O.D.I.,  Afolabi,  D.,  Damra,  F.M.,  On  the  Model  Stability  of  Perturbed  Cyclic 
Systems,  Advances  in  Control  and  Dynamics,  35, 137-164, 1990. 

4.  Afalobi,  D.,  Nwokah,  O.D.I.,  The  Frequency  Response  of  Mistuned  Periodic  Systems, 
Proc.  12th,  ASME  Biennial  Vibrations  Conference,  Montreal,  Canada,  September  1989. 

5.  Gilmore,  R.,  Catastrophe  Theory  for  Scientists  and  Engineers,  Wiley,  New  York,  1981. 

6.  Arnold,  V.I.,  Gusein-Zade,  S.M.,  Varchenko,  A.N.,  Singularities  of  Differentiable  Maps, 
Vol.  1,  Birkhauser,  Boston,  MA,  1985. 


-  16- 


7.  Arnold,  V.I.,  Catastrophe  Theory,  Springer  Verlag,  New  York,  1983. 

8.  Lu,  Y.C.,  Singularity  Theory  and  An  Introduction  to  Catastrophe  Theory,  Springer  Verlag, 
New  York,  1975. 

9.  Pierre,  C.,  Mode  Localization  and  Eigenvalue  Loci  Veering  Phenomena  in  Disordered 
Structures,  J.  Sound  and  Vibration,  126, 485-502k,  1988. 

10.  Pierre,  C.,  Cha,  P.D.,  Strong  Mode  Localization  in  Nearly  Periodic  Disordered  Structures, 
AIAA  Journal,  27,  227-241,  1989. 

11.  Pierre,  C.  Dowell,  E.H.,  Localization  of  Vibrations  by  Structural  Irregularity,  J.  Sound 
and  Vibration,  114,  549-564, 1987. 

12.  Afolabi,  D.,  Vibration  of  Mistuned  Bladed  Disk  Assemblies,  Ph.D.  Thesis,  University  of 
London,  England,  1983. 

13.  Davies,  P.J.,  Circulant  Matrices,  Wiley,  New  York,  1979. 

14.  Hamermesh,  M.,  Group  Theory  and  Its  Application  to  Physical  Problems,  Addison- 
Wesley,  London,  1962. 

15.  Arnold,  V.I.,  Singularity  Theory,  Cambridge  University  Press,  Cambridge,  1981. 

16.  Poston,  T.,  Stewan,  I.,  Catastrophe  Theory  and  Its  Applications,  Pitman,  London,  1978. 

STATEMENT  OF  WORK 

The  principal  aim  of  the  proposed  research  is  to  carry  out  an  in-depth  mathematical  and 
numerical  investigation  of  the  dynamics  of  mistuned  cyclic  systems,  by  use  of  some  new  and 
extremely  powerful  topological  theory  of  dynamics,  and  to  develop  simple  control  schemes  for 
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preventing  unacceptable  vibration  characteristics  in  such  systenas.  To  accomplish  this  task,  we 
will: 

(i)  Identify  the  topological  structure  of  nominally  mned  bladed-disk  assemblies,  the  order  of 
the  degeneracy  in  the  natural  frequencies,  the  minimum  number  of  canonical  parameters 
needed  to  unfold  the  degeneracy,  and  the  classification  of  the  bifurcation  set  in  the 
parameter  space. 

(ii)  Use  the  Jordan-Amold  canonical  structure  theory  to  completely  characterize  all  the  blade 
motion  forms  expected  when  a  given  nominally  tuned  system  is  generically  mistuned. 

(iii)  Relate  the  canonical  unfolding  parameters  to  the  disk  assembly  elements  of  mass, 
generalized  damping  and  generalized  stiffness;  and  hence  determine  which  mistuning 
parameters  or  combinations  thereof,  govern  the  escalation  of  forced  response  amplitudes 
and/or  unacceptable  blade  motions. 

(iv)  Employ  the  control  methodologies  of  either  entire  eigenstructure  assignment  or  quadratic 
optimization  to  deliberately  mistune  the  blade  assembly  passively  so  that  eigenvalue 
degeneracy  under  slight  parameter  variations  are  avoided  and  at  the  same  time  the 
parameter  combinations  which  lead  the  assembly  to  unacceptable  blade  motions  (the 
bifurcation  set)  are  never  allowed  to  occur. 

(v)  Carry  out  a  thorough  numerical  simulation  on  typical  nominal  and  peniu-bed  bladed-disk 
assemblies  to  verify  and  validate  the  predictions  of  the  new  topological  theory. 

The  above  will  set  the  stage  for  a  controlled  laboratory  hardware  experimental  verification, 
which  we  hope  to  undertake  in  a  follow-up  project. 
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ABSTRACT 

In  this  WOTk,  we  consider  the  dynamics  of  linear  mechanical  systems  possessing  geometri¬ 
cal  symmetry  subject  to  differential  or  small  parameter  variations.  The  machinery  of  group 
theory  including  the  irreducible  group  representations,  and  the  consideration  of  representations 
to  which  the  translational,  rotational  and  vibrational  modes  belong,  allow  us  to  predict  apriori, 
the  number  and  the  order  of  degenerate  eigenvalues  in  the  symmetric  system.  By  considering 
the  resultant  Hamiltonians  of  the  perturbed  symmetric  system,  we  show  further  the  effects  of  the 
perturbadons  on  the  eigenvalues  and  their  degeneracies.  Since  the  vibradon  modes  of  systems 
with  degenerate  eigenvalues  are  known  to  display  sensidve  dependence  on  parameters,  we  may 
use  these  techniques  to  identify  in  principle  the  possibility  of  maximum  vibradon  amplitudes  and 
where  they  are  likely  to  occur.  Applications  of  these  ideas  include  the  mistuned  turbine  rotor 


bladed  disk  assemblies. 
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LIST  OF  SPECIAL  SYMBOLS 

r  (R)  Reduced  cartesian  representation  of  a  group  element  R. 
r  ^  (R)  Reduced  translational  representation  of  a  group  element  R. 

Reduced  Rotational  representation, 
r Reduced  Vibrational  representation. 
r(R),  r^R)  Matrix  representations. 
r'(R)  i*  representation. 

riiv(R)  Matrix  element  of  the  jx*  row  and  the  v*  column  of  the  matrix  representing  the 

group  element  R  in  the  i*  representation. 

r^v  (R)  Complex  conjugate  of  r{iv  (R). 

X‘(R)  Character  of  a  group  element  R  in  the  i*  matrix  representation, 

a;  Number  of  times  F  ‘(R)  appears  in  the  reducible  representation, 

e  Row  vector. 

M  Mass  matrix. 

K  Stiffness  matrix. 
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1.  INTRODUCTION 

Eigenvalues  and  eigenvectors  of  a  vibrating  system  are  important  for  characterizing  its 
dynamical  response.  The  eigenvalues  are  related  to  natural  frequencies  whereas  the  eigenvectors 
correspond  to  special  forms  of  displacements  when  vibrating  at  a  natural  frequency.  Exact 
evaluation  of  the  eigenvalues  of  higher  order  vibrating  systems  in  general  involves  considerable 
effort  and  is  time  consuming.  Most  cyclic  symmetric  systems  possess  degenerate  eigenvalues 
[1,2].  Systems  with  degenerate  eigenvalues  arc  expected  to  display  severe  sensitive  dependence 
on  parameters  [3]  that  destroy  the  symmetry  or  degeneracy. 

In  the  eigenvalue  problem  if  there  is  any  symmetry  of  the  system,  the  application  of  group 
theory  enables  us  to  decide,  at  the  outset,  exactly  the  number  of  distinct  eigenvalues  together 
with  their  respective  degrees  and  degeneracies. 

By  considering  the  symmetry  operations  of  the  physical  system  at  the  equilibrium  points, 
the  representing  group  can  be  formulated.  Using  group  theoretical  ideas,  we  can  predict  apriori 
the  degeneracy  of  the  eigenvalues.  This  is  accomplished  by  the  use  of  the  irreducible  represen¬ 
tations  of  this  group  which  is  obtained  by  using  the  orthogonality  theorem  and  the  reduction  for¬ 
mula  [3].  Once  the  irreducible  representations  are  known,  we  can  find  the  translational,  rota¬ 
tional  and  vibrational  modes  of  the  system.  These  results  are  well  known  in  the  physics  litera¬ 
ture  on  group  theory  but  have  not  been  used  sufficiently  effectively  in  the  vibration  community. 
The  essential  purpose  of  this  work  is  to  summarize  some  of  these  results  and  show  some  applica¬ 
tions  as  they  relate  to  the  symmetric  bladed  disk  assemblies. 

In  general,  a  reduced  canesian  representation  of  a  group  element  R,  r‘^(R),  can  be  writ¬ 


ten  as 
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r®-"  (R)  =r‘^  (R)®r'-*  (R)®r''-‘  (R), 

where  r‘^(R),  r*^(R),  and  r''*^(R)  are  translational,  rotational  and  vibrational  reduced  carte¬ 
sian  representations  of  a  group  element  R.  In  three  dimensional  vibrating  systems,  six  of  the 
normal  modes  belong  to  the  zero  frequency  modes  and  correspond  to  pure  translations  and  pure 
rotations.  Since  we  are  primarily  interested  in  vibrational  modes  (non  zero  frequency)  the  zero 
frequency  modes  are  not  discussed  fmther  in  this  work. 

Some  of  the  symmetry  of  the  physical  system  may  be  lost  once  the  system  is  subjected  to 
parameter  perturbations.  The  perturbed  system  may,  however,  still  possess  some  symmetry 
which  may  be  considered  to  be  a  subgroup  of  the  group  characterizing  the  unperturbed  or  origi¬ 
nal  system.  Applying  group  theoretic  ideas  now  to  this  subgroup  we  can  predict  the  splitting  of 
the  degeneracy  of  eigenvalues.  As  a  result  we  can  see  whether  the  degeheracy  of  some  of  the 
eigenvalues  has  or  has  not  been  removed. 

Ideas  similar  to  the  ones  proposed  in  this  woric  were  used  by  R.  Perrin  [4]  in  1971  for  a 
thin  circular  ring.  In  his  paper  group  theoretical  arguments  were  applied  to  a  ring  where  pertur¬ 
bation  was  applied  in  the  form  of  equal  masses  attached  to  the  ring  at  the  vertices  of  an  inscribed 
regular  n*  order  polygon.  Further,  eigenfrequencies  and  eigenfunctions  for  the  unperturbed  ring 
were  assumed  to  be  known  apriori.  Knowledge  of  these  degenerate  pairs  of  eigenfunctions  was 
used  to  find  the  characters  of  each  irreducible  representation  of  the  corresponding  Dn  group  for 
the  perturbed  system.  In  the  present  work,  group  theoretical  techniques  are  developed  without 
apriori  assuming  any  knowledge  of  either  the  eigenfrequencies  or  the  eigenfunctions  for  the  cir- 
culant  symmetric  system.  Also,  parametrically  perturbed  cases  were  not  discussed  by  Perrin. 
Parametric  perturbations  arc  important  in  turbine  blade  vibration  problems  where  a  slight  pertur¬ 
bation  can  lead  to  loss  of  cyclic  symmetry,  which  in  turn  can  induce  rogue  blade  failure  under 
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cenain  circumstances  [5,6]. 

2.  THE  GROUP  THEORETICAL  CONCEPTS 

We  first  define  some  standard  terms  from  the  literature  on  group  theory.  One  of  the  stan¬ 
dard  references  is  the  text  by  Hamermesh  [1].  Following  definitions  and  theorems  are  obtained 
from  [2]. 

2.1  Definition  1.  Symmetry  operations:  All  the  operations  which  leave  a  system 
configuration  unchanged  are  called  symmetry  operations. 

In  physical  terms,  this  refers  to  the  movement  of  a  system  in  such  a 
way  that  it  interchanges  the  positions  of  various  particles  of  the  system  but 
results  in  the  system  looking  exactly  the  same  as  before  the  symmetry 
operation.  For  instance  some  of  the  symmetry  operations  are  defined  as 
follows; 

E:  Identity.  The  system  is  not  rotated  at  ail  or  rotated  by  27C  about  any 

axis. 

Cn  Rotation:  This  is  an  operation  which  effects  rotation  through  an 
angle  27c/n  about  an  axis,  fixed  in  space,  where  n  is  an  integer.  In 
addition  we  can  have  which  is  €„  raised  to  the  power  k,  that  is, 
a  rotation  through  an  angle  27tk/n  about  the  same  axis.  CH  is  a  rota¬ 
tion  through  an  angle  2n  and  is  the  identity  operation,  since  a  rota¬ 
tion  through  2n  leaves  the  object  unchanged,  n  is  known  as  the 
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Definidon  2. 


Definition  3. 


multiplicity  of  the  axis,  and  the  latter  is  called  on  n  fold  axis.  If 
n=!2,3,...  then,  respectively,  we  get  2-fold,  3-fold...  axes.  If  a  sys¬ 
tem  has  more  than  one  axis  of  symmetry  then  the  axis  with  the 
highest  value  of  n  is  called  the  principal  axis. 

Group:  A  set  of  elements  (a,b,c..}  is  called  a  group  G,  if  a  multiplication 
rule  is  defined  for  any  two  elements  so  that  the  product  ab  has  a  definite 
meaning  and  the  following  four  postulates  are  satisfied: 

1.  Closure:  If  a  and  b  belong  to  the  set,  then  ab  also  belongs  to  the  set. 

2.  Associativity:  a(bc)  =  (ab)  c. 

3.  There  exists  the  identity  element  e  such  that  ae  =  ea  =  a  for  any  a 
belonging  to  G. 

4.  There  exists  the  inverse  element,  i.e.,  for  each  element  a,  there  is  a 
corresponding  element  b  such  that  ab  =  ba  =  e.  b  is  called  the 
inverse  element  of  a  and  is  denoted  by  b  =  a~^ . 

Dn  group:  This  group  concerns  a  system  possessing  one  n-fold  axis  called 
the  principal  axis  and  n  2-fold  axes  symmetrically  placed  in  a  plane  per¬ 
pendicular  to  the  principal  axis.  The  n-fold  axis  provides  the  n  elements 
of  the  cyclic  group  Q,.  The  group  also  contains  one  C2  element  provided 
by  every  perpendicular  2-foId  axis  where  we  do  not  count  C2  =  E  because 
E  only  occurs  once  in  a  set  of  group  elements  and  it  has  already  appeared 
in  the  n  elements  of  Q,.  The  group  Dn  therefore  contains  a  total  of  2n  ele¬ 


ments. 
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Definidon  4.  Equivalent  and  reducible  representations:  Two  representadons  are  said 

to  be  equivalent  if  the  two  matrices  representing  any  element  R  of  the 
group  are  related  by  the  equadon 

r(R)T,  (1) 


where  T  is  any  nonsingular  square  matrix  (operator).  However,  if  there 
does  not  exist  any  matrix  T  which  transforms  P '  (R)  into  F  (R),  then 
r '  (R)  and  F  (R)  are  said  to  be  inequivalent. 


A  reduced  representadon  of  a  group  element  R,  F^*^  (R),  is  com¬ 
posed  of  two  or  more  irreducible  representadons: 


F«^(R)  = 


F‘  (R)  0 

0  Pj  (R) 


(2) 


We  write  this  by  using  the  symbol  ©  and 


F^«*  (R)  =  F‘  (R)©rJ  (R) 


(3) 


Definidon  5.  X*(*’)*  Character  of  a  group:  The  character  of  a  group  element  R  in  the  i* 

matrix  representadon  of  the  group  element  is  the  trace  (sum  of  diagonal 
terms)  of  the  matrix. 


2.2.1  Orthogonality  Theorem:  All  the  vectors  formed  by  the  inequivalent  irreducible  unitary 
representadons  are  orthogonal  to  each  other,  or: 

I  rlTv  (R)  Fjiv  (R)  =  f  5ij5„,'5w'. 

R  *' 


where  i  and  j  denote  the  representation,  p.  and  \i'  denote  rows  of  the  matrix  elements  and  v  and 
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V '  denote  the  columns  of  the  matrix  elements,  g  is  the  order  of  the  group  and  is  the  dimen¬ 
sionality  of  the  i***  representation.  r|iv  (R)  is  the  matrix  element  of  the  row  and  column 
of  the  matrix  representing  the  group  element  R  in  the  i***  representadon,  and denote  its  complex 
conjugate  equivalent  to  F  ‘  (R) 

2.2.2  Character  table:  The  character  table  is  formed  by  considering  the  characters  of  group 
elements.  The  character  of  a  group  element  is  important  because  the  character  is  unaltered  by  a 
similarity  transformation.  On  the  other  hand  since  the  character  of  equivalent  irreducible 
representations  are  identical  a  table  of  characters  is  a  unique  way  to  characterize  a  group.  The 
general  form  of  a  character  table  is: 


where 

c(r);  the  nature  of  elements  in  the  class, 

Ni :  number  of  elements  in  the  class,  and 
r‘:  i'*’ irreducible  representation. 

The  character  tables  of  the  groups  can  be  obtained  by  the  application  of  the  following  rules: 
( 1 )  Number  of  inequivalent  representations  is  equal  to  the  number  of  classes. 
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(2)  £  if  =  g>  n-  number  of  irreducible  representations. 

1=1 

(3)  2  iC*  (R)  JC*  (R)  =  g  where  x‘(R)  denote  the  character  of  a  group  element  R  in  the  i* 

R 

matrix  representations. 

(4)  ix''(Ci,)x‘(Ck)=f^Sb.- 

In  labeling  the  rows  of  the  character  table,  the  following  standard  notation  is  used. 

(1)  One-dimensional  representations  are  labeled  as  A  if  the  character  of  the  elements  about 
the  principal  rotation  axis  are  +1  for  all  k,  and  as  B  if  the  characters  Cn  are  (-1)*^  for  all  k. 

(2)  If  a  group  has  more  than  one  A  or  B  representation  they  are  given  subscripts  1  and  2 
according  to  whether  the  character  is  +1  or  -1  in  the  column  representing  a  rotation  or 
improper  rotation  about  an  axis  other  than  the  Piincipal  axis.  For  example,  in  the  groups 
Dn  a  representation  is  given  a  subscript  1  if  the  character  under  C2  about  the  axis  is  +1  and 
2  if  it  is-1. 

(3)  Two  dimensional  representations  are  labeled  E. 

(4)  Three  dimensional  representations  are  usually  labeled  T. 

Finally,  the  reducible  representations  are  used  to  obtain  the  irreducible  representations  by  the 
applications  of  the  formula 

aj  =  -INjXJ'(R)x'“*(R).  (4) 

8  R 
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3.  VIBRATIONS  OF  CYCLIC  MECHANICAL  SYSTEMS 

Many  vibrating  systems  possess  sufficient  symmetry  to  allow  us  to  use  group  theory  which 
reduces  the  amount  of  work  involved  in  the  calculations  and  also  furnishes  us  with  an  insight 
into  the  nature  of  the  vibrations.  Consider  n  nominally  identical  masses  connected  via  ground 
springs  kt  and  coupling  springs  kc  at  the  edges  of  an  inscribed  n*  order  regular  polygon  in  a  cir¬ 
cle  of  radius  r.  Rotational  symmetry  about  one  n-fold  axis  perpendicular  to  the  plane  of  motion 
and  n,  2-fold  axes  give  2n  number  of  elements  for  the  corresponding  symmetry  group.  In  fact, 

this  is  the  dihedral  group  Dn.  Dn  has  (2  +  --■■-■)  conjugate  classes  when  n  is  odd  and  (3  +  y) 

classes  when  n  is  even.  Utilizing  group  character  table  construction  rules  [1,2],  it  can  be  shown 
that  the  number  of  possible  degenerate  eigen  levels  are: 


n  even 

single  degenerate  levels  s  4 

i)  n  ^ 

double  degenerate  levels  =—  -  1  ’ 

n  odd 

single  degenerate  levels  =  2 
double  degenerate  levels  = 


(5) 


(6) 


Results  (i)  and  (ii)  imply  that  a  cyclic  symmetric  system  of  this  type,  at  worst,  can  have  double 
degenerate  vibrational  modes  for  any  finite  n. 

Once  the  system  is  subjected  to  a  random  parameter  perturbation,  some  of  the  symmetry 
may  be  lost  and  consequently,  we  get  a  new  group  which  is  likely  to  be  a  subgroup  of  the  ohgi- 
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nal  group.  Equations  (5)  and  (6)  can  be  used  to  get  a  qualitative  information  about  the  new 
degenerate  eigenvalues  of  the  perturbed  system.  Since  any  higher  order  cyclic  symmetric  sys¬ 
tem  of  this  type  will  have  at  most  doubly  degenerate  eigenvalues  it  is  sufficient  to  consider  an 
example  with  n  =  3. 

Consider  three  normally  identical  masses  connected  via  both,  ground  springs  kt  and  cou¬ 
pling  springs  kc,  at  the  edges  of  an  inscribed  isosceles  triangle,  in  a  circle  of  radius  r.  We  wish 
to  study  the  number  and  degeneracy  of  the  eigenvalues  of  such  systems  under  random  differen¬ 
tial  perturbations  in  the  elements  {kf  k^.,  m). 

As  defined  in  Section  2,  a  symmetry  operation  is  one  which  leaves  the  undistorted  system 
indistinguishable  from  its  previous  orientation.  Such  an  operation  interchanges  equivalent 
masses.  However,  in  the  vibrational  state,  the  system  is  in  a  distorted  configuration  and,  when 
the  symmetry  operation  is  performed  on  the  distorted  mass  the  effect  is  the  same  as  that  obtained 
by  interchanging  displacement  vectors  amongst  equivalent  masses.  Therefore  we  can  define  the 
action  of  a  symmetry  operation  for  each  mass  in  a  distorted  system  to  be  a  displacement  through 
vector  Xj  from  its  equilibrium.  When  a  symmetry  operation  is  applied  we  can  assume  that  the 
mass  positions  remain  invariant.  Furthermore  the  symmetry  operation  can  have  no  effect  on  the 
potential  or  kinetic  energy  of  the  system,  or  even  the  angles  between  the  connections.  Conse¬ 
quently  the  quadratic  forms  of  the  kinetic  energy  T  and  the  potential  energy  V  remain  invariant 
under  the  action  of  the  group  transformations.  Group  theory  can  thus  be  used  to  determine  and 
classify  the  normal  modes  of  the  vibrating  system. 

We  begin  with  the  3N  dimensional  representation  of  the  group  of  symmetry  operations  of 
the  undistorted  system.  By  reducing  this  representation  using  the  character  table  for  the 
corresponding  group,  (and  reduction  of  reducible  representations)  we  can  determine  the 
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intducible  representations  to  which  the  3N  translational,  rotational  and  vibrational  modes 
belong.  Also  we  can  immediately  find  the  degeneracy  of  each  normal  mode.  In  addition,  by 
considering  the  symmetrized  basis  S,  we  can  bring  the  mass  matrix  M  and  the  stiffness  matrix  K 
into  block  form  and  thus  greatly  simplifying  the  solution  of  the  characteristic  or  the  frequency 
equation. 

We  apply  these  group  theoretic  techniques  to  the  problem  shown  in  Figiue  1  [4].  This  sys¬ 
tem  belongs  to  the  group  D3  which  contains  the  symmetry  operations  E,  C^,  C^,  C|,  C^,  and 
Cf.  Now  by  applying  these  operations  to  the  nine  canesian  coordinates  xi,  X2,  X3,  y^,  y2,  y3,  zj, 
Z2,  Z3  we  obtain  their  nine  dimensional  reducible  representations.  This  can  be  accomplished  by 
finding  the  corresponding  matrix  representation,  and  using  the  equation 

X'*r(R)X,  (7) 

where  R  is  the  appropriate  symmetry  operation,  and  X  and  X'  denote  the  nine  dimensional  vec¬ 
tors  representing  the  cartesian  coordinates  in  original  and  transformed  planes  respectively.  For 
example,  under  the  operation  of  the  system  configuration  in  Figure  1  is  transformed  to  the 
configuration  in  Figure  2.  This  also  clearly  indicates  the  manner  in  which  the  coordinates 
undergo  rotation.  In  Figure  2,  the  axes  zj,  Z2  and  Z3  are  pointing  out  of  the  plane  of  the  paper. 
The  new  and  the  old  coordinates  are  related  by  the  relations 

(xj ,  yi  .Zj )  -» (x'l ,  /i ,  z'l )  ■  (-  X2  sin  30  -  y2  cos  30,  X2  cos  30  -  y2  sin  30,  Z2) ,  (8) 

(X2,  y2.  Z2)  -*  (x'2,  z'2)  ■  (-  X3  sin  30  -  y3  cos  30,  X3  cos  30  -  y2  sin  30,  Z3) ,  (9) 

(^3.  y3,  Z3)  (x'3,  /3,  z'3)  ■  ( -  X3  sin  30 - y3  cos  30,  X3  cos  30 - y2  sin  30,  Z3) .  (10) 

In  matrix  form,  equations  (8)  -  (10)  are  represented  as 
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(11) 


where  the  coefficient  matrix  is  the  matrix  representation  of  the  group  element  Cj.  The  character 
(trace)  of  the  matrix  is  then 

X"(Cl)  =  0.  (12) 

Since  and  Cl  belong  to  the  same  class,  X®  (C^)  =  0.  It  is  clear  that  x®  (E)  =  9.  In  a  similar 


manner 

x'(C5)=-i, 

X'(C!)  =  -1,  (13) 

x‘(C!)=-i. 

Now  using  the  character  table  for  D3  (Table  1)  and  equation  (4),  we  can  determine  the 
irreducible  components  off®  (R)  as: 


Table  1 
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Here,  Ai ,  A2  represent  one  dimensional  representations,  E  represents  two  dimensional 
representations,  and  Tx,  Ty,  Tj  and  Rx,  Ry,  Rz  represent  unit  translational  and  rotational  vectors 
respectively.  Consequently 

=T(^^i+2(0xl)-t-3(lx-l))=l. 

o 

aA,  =-|(9x  l  +  2(0x  l)  +  3(-lx-l))  =  2, 
aA,  =  -^  (9  X  2  +  2  (0  X  -1)  +  3  (1  X  0))  =  3  . 

Hence  (R)  =  Aj  G  2A2  G  3E.  Since  the  system  has  three  masses,  N  =  3  and  there  are  9 
degrees-of-freedom  for  the  system.  Since  3N  -  6  =  3,  we  have  only  three  vibrational  modes,  the 
other  six  corresponding  to  zero  frequency  modes  and  to  pure  rotations  and  pure  translations. 
This  can  be  seen  from  the  character  table, 

r‘(R)  =  A2GE, 
r  (R)  =  A2GE. 

Since  we  are  primarily  interested  in  the  vibrational  modes  of  the  system  and 

r®-  (R)  =  Ai  G2A2G3E, 
r‘(R)Gr(R)  =  2A2G2E, 

we  get  r'^"'  (R)  =  Ai  G  E. 

Since  A\  is  one  dimensional,  the  vibrating  system  has  one  non  degenerate  eigenvalue  and 
since  E  is  2  dimensional  there  is  a  one  degenerate  eigenvalue.  These  results  are  consistent  with 
the  exact  eigenvalues  given  in  appendix  A. 


Using  group  theory  arguments  to  predict  the  number  of  degenerate  eigenvalues  becomes 
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very  useful  as  the  order  of  the  system  increases  since  the  exact  calculation  of  eigenvalues  for 
such  high  order  systems  becomes  increasingly  burdensome. 

4.  THE  PERTURBED  SYSTEM 

Suppose  that  the  Hamiltonian  of  the  unperturbed  system  is  Hq.  Then  Hq  is  invariant  under 
its  symmetry  group  G.  Suppose  funher  that  the  system  is  subjected  to  a  perturbation  with  Ham¬ 
iltonian  V.  The  perturbed  Hamiltonian  H  =  Hq  +  V,  will  then  have  a  symmetry  group  which  is 
necessarily  a  subgroup  of  G.  Two  possible  cases  arise. 

CASE  I 

If  the  perturbation  V  has  symmetry  at  least  as  great  as  Hq.  the  group  G  will  still  be  the 
symmetry  group  of  the  total  Hamiltonian  H.  In  this  case  the  possible  types  of  eigenvalues  will 
be  unchanged  by  the  perturbation.  In  fact  no  splitting  of  degenerate  levels  occurs. 

CASEn 

If  the  perturbation  V  has  symmetry  lower  than  Hq,  the  total  Hamiltonian  H  will  have  a 
symmetry  group  G*  which  is  a  subgroup  of  G.  This  subgroup  G'  is  invariant  under  the  perturba¬ 
tion.  Because  of  the  perturbation,  some  of  the  degenerate  eigenvalues  may  split.  This  can  be 
explained  by  using  group  representation  theory. 

For  a  given  representation  D(G)  of  the  group  G,  we  now  obtain  the  invariant  subgroup  G'  . 
Even  if  D(G)  is  an  irreducible  representation  of  G.  the  representations  of  G^  which  we  derive  in 
this  way  may  be  reducible.  In  other  words,  even  though  we  cannot  find  a  subset  of  the  basis  vec¬ 
tors  of  D(G)  which  is  invariant  under  all  transformations  of  the  group  G,  we  may  be  able  to  find 
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a  subspace  which  is  invariant  under  all  transformations  belonging  to  the  eigenvalue  X  from  a 
basis  of  an  irreducible  representation  of  G.  This  representation  may  be  reducible  for  the  sub¬ 
group  G'  .  The  perturbation  V  will  then  split  the  level. 

We  now  apply  the  above  ideas  and  show  the  appropriate  methodology  in  the  context  of  the 
three  mass  system. 

4.1  Ground  Spring  Perturbed  System 

First  we  consider  the  situation  when  one  of  the  ground  springs  is  perturbed.  The  unper¬ 
turbed  system  can  be  represented  by  the  group  D3  =  {E,  C3,  C^,  C*.  Cb,  Cg).  Once  the  system 
is  subjected  to  a  ground  spring  perturbation,  for  this  particular  system  G'  =  {E,  C,}  is  the  invari¬ 
ant  subgroup. 

The  character  table  for  G^  =  {E,  C, ),  is  as  follows: 


Considering  the  part  of  the  character  table  of  D3  which  refers  to  the  operations  of  the  sub¬ 
group  =  {E,  C, },  we  have 


Utilizing  (4),  the  irreducible  components  are  then  given  by 
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^  (1  X  2  X  1  +  1  X  0  X  1)  =  1 , 

A 

=  y  (1  X  2  X  1  +  1  X  0  -1)  =  1 . 

Thus,  the  doubly  degenerate  E  level  of  the  unperturbed  system  splits  into  single  levels  A'  and  A" 
of  under  the  ground  spring  perturbation.  As  a  result,  degenerate  eigenvalue  of  the  perturbed 
system  separates.  Hence,  for  this  particular  system  we  get  three  distinct  eigenvalues.  Coupling 
spring  perturbation  leads  to  a  case  where  there  is  no  invariant  subgroup  left  and  consequently, 
group  theoretical  arguments  do  not  work  for  this  particular  situation.  We  conjecture  that  this  is 
indicative  of  those  cases  where  perturbations  do  not  lead  to  radical  changes  in  the  eigenvector 
directions. 

4.2  The  Mass  Perturbed  System 

A  system  consisting  of  three  particles,  two  with  mass  m  and  the  other  with  mass  M,  is 
illustrated  in  Figure  3,  Considering  rotational  symmetry  of  the  system,  we  can  see  from  Figure  3 
that  this  system  belongs  to  the  group  D2 

D2  =  {E.  C2). 

The  character  of  each  of  the  elements  of  the  group  can  be  obtained  from  the  reducible  represen¬ 
tation  whose  matrix  representation  is  obtained  by  the  use  of  the  coordinate  transformation 

=  r  (R)  X. 

Performing  these  operations,  and  following  the  steps  along  the  lines  of  work  in  section  3,  we  can 
show  that  the  reducible  characters  of  each  element  are 


X‘^(C2)  =  -l, 
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X'(E)  =  9. 

Using  the  character  table  of  D2  and  (4),  the  irreducible  components  of  the  reducible  representa¬ 
tions  can  be  then  determined  as 


D2 

E 

C2 

A 

1 

1 

z 

Rz 

B 

1 

-1 

x,y 

Rx»  Ry 

aA  =  y  [9xl-lxlj=4, 
aa  =  y  [9X  1+-1  X  -  ij  =  5. 


Therefore, 

(R)  =  4A@5B.  (14) 

In  general,  T  (R)  =  T  (R)  ®  (R)  ©  P’'"'  (R). 

By  placing  a  coordinate  system  XYZ  at  the  center  of  mass,  translations  and  rotations  can  be 
represented  as  shown  in  Figure  4.  The  representations  for  C2  and  E  are  given  by 

X‘(C2)  =  -1  .  X‘(E)  =  3. 

X'(C2)=-l  .  X^(E)  =  3. 

By  the  application  of  the  aj  equations  (4),  we  get 


(i)  Translation 
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ak  =  Y[lx3-lxl]  =  l, 
ak  =  4-llx3-lx-l]  =  2. 

(ii)  Rotation 

ak  =  j[lx3-lxl]  =  l, 
afe=Y[lx3-lx-l]  =  2. 

Therefore, 


r‘~‘  (R)  =  A@2B  , 

(15) 

r"**  (R)  =  A  @  2  B  . 

(16) 

Hence  from  equations  (14)-(16),  we  get  the  result  that 

r''“‘ (R)  =  2A@B  .  (17) 

This  shows  that  there  are  three  vibrational  modes  and  each  eigen  level  is  non-degenerate  since  A 
and  B  are  one  dimensional  representations. 

5.  SUMMARY  AND  CONCLUSIONS 

This  work  uses  results  from  group  theory  and  applies  it  to  perturbed  cyclic  symmetric 
vibratory  systems.  It  is  shown  that: 

a.  The  number  and  order  of  degenerate  eigenvalues  in  a  symmetric  system  can  be  predicted 
apriori  by  using  group  theory  without  explicitly  determining  the  eigenvalues. 

b.  Cyclic  symmetric  vibrating  systems  possess  degenerate  eigenvalues  for  n  >  2.  For  strong 


-21  - 


coupling,  these  eigenvalues  occur  in  doubly  degenerate  pairs  and  in  single  nondegenerate 
levels. 

c.  Random  parameter  perturbations  may  partially  or  totally  destroy  the  symmetry  of  the  sys¬ 
tem.  Accordingly  these  perturbations  lift  some  of  the  degeneracy  of  eigenvalues.  As  a 
result,  eigenvalue  loci  veering  [7]  occurs  when  the  parameters  are  continuously  varied.  This 
may  also  lead  to  a  mode  localization  or  rapid  variation  in  the  eigenfunctions  [7]. 
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Appendix  A 

Equations  of  motion  for  the  system  in  Figure  2  can  be  written  in  the  form: 


Mx  +  Kx  =  0  ,  where 


o 

o 

S 

a  -kc  -kc 

M  = 

0  m  0 

.  K= 

-kc  a  -kc 

0  0m 

-kc  -kc  a 

a  =  2kc  +  . 


The  eigenvalues  are  determined  by 


det  1 K  -  0)^  M  1  =  0 . 


Therefore, 


-kt  3kc  3kc 

(of  =  — r-  + -  ,  ©2  =  — T  + - » 

mr  in  mr^  m 


This  shows  that  the  cyclic  symmetric  system  has  a  double  degenerate  eigenvalue  and  a  sin¬ 
gle  nondegenerate  eigenvalue. 


LIST  OF  FIGURES 

Figure  1.  Model  representing  the  unperturbed  three  bladed  disk  assembly  in  its  identity 
orientation. 

Figure  2.  Model  of  figure  1  counter/clockwise  rotated  by  2jt/3  radians  about  the  center  of  the 
system,  showing  change  in  cartesian  coordinate  axes. 

Figure  3.  Mass  perturbed  system  in  its  identity  orientation. 

Figure  4.  The  representation  of  the  three  translations  and  the  three  rotations  in  an  XYZ 
coordinate  system. 


kc 


Figure  2.  Model  of  figure  1  counter/clockwise  rotated  by  2)1/3  radians  about  the  center  of  the 
system,  showing  change  in  cartesian  coordinate  axes. 
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In  recent  years  there  has  been  tremendous  interest  in  the  vibrations  and  structural  dynamics 
community  in  the  phenomenon  of  mode  localization.  This  interest  stems  from  the  recognition 
that  large  systems  composed  of  nominally  identical  subsystems  inevitably  involve  minor  devia¬ 
tions  from  the  idealized  structures  and  these  disorders  or  perturbations  can,  under  appropriate 
conditions,  cause  disproportionately  large  deviations  from  the  predicted  behavior  in  the  nominal 
or  idealized  system  modes.  Imponant  technical  applications  of  these  include  mistuned  bladed 
disk  assemblies  [1]  and  large  space  structures  (21. 

It  is  well  understood  by  now  that  the  presence  of  small  irregularities  in  nearly  periodic  struc¬ 
tures  may  inhibit  the  propagation  of  vibration  and  localize  the  vibration  modes.  Depending  on 
the  magnitude  of  perturbations  (disorder)  and  on  the  strength  of  internal  coupling  between  the 
subsystems,  the  mode  shapes  may  undergo  dramatic  changes  to  become  strongly  localized  when 
small  perturbations  are  introduced,  thereby  confining  the  energy  associated  with  a  given  mode  to 
a  small  geometric  region.  This  phenomenon  is  referred  to  as  mode  localization.  Pierre  [3] 
showed  that  strong  mode  localization  and  eigenvalue  curve  veering,  are  two  manifestations  of 
the  same  phenomenon.  Therefore,  the  investigation  of  the  curves  of  the  eigenvalues  or  natural 
frequencies  in  the  neighborhood  of  the  ordered  state  is  sufficient  for  determining  the  occurrence 
of  strong  mode  localization.  That  eigenloci  veering  phenomenon  can  occur  in  disordered  struc¬ 
tures  under  certain  conditions  has  been  e.xplained  qualitatively  using  geometric  arguments  m  [4]; 
where  it  was  also  hinted  that  quantitative  results  should  be  obtainable  by  the  use  of  singular  per¬ 
turbation  analysis. 

For  mistuned  linear  chains,  Pierre  (3|  and,  Pierre  and  Dowell  [5]  showed  that  the  straight  for¬ 
ward  expansion  in  terms  of  mistuning  parameters  breakdown  in  the  case  of  weak  coupling.  This 
arises  because  the  idealized  system  that  is  being  penurbed  has  natural  frequencies  with 
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multiplicity  >  1.  They  then  developed  a  so-called  "modified"  penurbation  technique  which  pro¬ 
vided  a  good  approximation  to  the  exact  eigenfrequencies  and  showed  good  agreement  with 
experimental  results. 

In  the  present  note  we  show  that  the  singularity  causing  the  breakdown  of  the  straight  for¬ 
ward  expansion  can  be  analyzed  by  the  uell  developed  singular  perturbation  techniques  [6|  and 
an  appropriate  asymptotic  expanstiin  for  the  eigenfrequencies  can  be  constructed  which  provides 
a  correct  qualitative  and  good  quantitative  approximation.  In  order  to  explain  the  ideas  and  to 
keep  the  algebraic  manipulations  to  a  minimum,  the  attention  is  focused  on  the  now  standard 
example  [3]  of  the  coupled  penduli  shown  in  Figure  1. 

The  basic  idea  of  the  technique  is  the  following:  by  applying  the  regular  perturbation  tech¬ 
nique  to  the  characteristic  equation  F(X,  e,  6)  =  0,  of  the  system,  we  can  obtain  algebraic  expres¬ 
sions  for  the  natural  frequencies  (eigenvalues)  as  a  power  series  in  the  small  parameter  or  penur¬ 
bation  (say  5).  The  coefficients  of  the  power  series  are  dependent  on  the  second  parameter  e  and 
these  expansions  are  valid  for  all  values  of  e  so  long  as  no  singularities  arise.  Singularities  occur 
for  values  of  e  where  the  eigenfrequencies  lose  their  smoothness  and  it  is  said  that  the  expansion 
is  not  uniformly  valid  for  all  e.  Away  from  the  singular  parameter  (e)  values,  the  straightforward 
expansions  are  good  approximations  and  are  called  the  "outer  expansions The  neighborhood 
of  the  singular  parameter  point  i>  then  stretched  or  rescaled  in  terms  of  a  new  parameter  so  as  to 
remove  the  singularity.  The  expansions  in  terms  of  the  new  parameter  is  valid  only  in  the  neigh¬ 
borhood  of  the  singular  point  and  is  c.illed  the  "inner  expansion".  The  inner  and  the  outer  solu¬ 
tions  can  be  matched  where  their  domains  of  validity  overlap  and  then  a  composite  expansion 
can  be  constructed  which  is  valid  uniformly  throughout  the  function  domain  for  all  values  of  the 


parameter  e. 
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Considcr  the  system  of  two  weakly  coupled  penduli  system  as  shown  in  Figure  1.  The  two 
imponant  parameters  are  the  dimensionless  coupling  between  pendulums  =  (k/m)/(g/l),  and 
the  dimensionless  length  change  The  corresponding  eigenvalue  problem  generated  by  the 
above  system  is  given  by; 


1+R^  -R^ 

-R^  R"+(l+^/r' 


<t)  =  xo 


ll.li 


where  R^  =  cOk/cOg  .  cok  =  k/'m  ,  co"  =  g//  . 


This  eigenvalue  problem  results  in  the  following  characteristic  equation: 


F(>.,e,5)  =  X^-i  i+2&4-^)?l  +  —  +  5  =  0. 

l+e  1+e 


(1.2) 


where  A/  =  e  and  R^  =  5. 


We  can  express  the  solutions  to  (1.2)  as  regular  functions  of  the  parameters  A/,  6  as  follows: 


X.1  (A/)  =  1  +  5  + 


5^  +0(6^)  =  Xf  (A/), 


(1.3) 


X2(A/)  = 


I 


I  +A/ 


+  5- 


5^ +0(5^)  =  XI  (A/) 


(1.4) 


The  expressions  (1.3),  (1.4)  are  the  regular  expansions  of  the  eigenvalue  problem  for  small  cou¬ 
pling.  When  A/  0,  Xi  and  Xi  become  unbounded  and  the  continuity  of  the  eigenvalues  with 
respect  to  the  perturbation  A/  breaks  down,  as  shown  in  Figure  2.  Each  eigencurve  has  two 
branches,  one  valid  for  A/  >  0  and  the  other  for  A/  <  0.  These  branches  are  indicated  by  the 
superscripts  ’+’  and  which  correspond  to  A/  >  0  and  A/  <  0,  respectively.  Note  that,  since  A/ 
and  R^  =  6  have  been  treated  as  two  independent  parameters,  we  have  no  control  over  expres- 
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sions  (1.3)  and  (1.4)  in  the  limiting  process  when  A/  -4  0  and  5  — >  0.  B'-  forming  the  "inner 
e.xpansion"  however,  we  can  find  an  e.xact  relation  between  these  two  parameters  by  talking  into 
consideration  the  nature  of  the  singularity.  Then  A/  and  5  become  dependent  parameters. 
.Asymptotically  matching  the  inner  and  the  outer  expansions,  then  gives  the  composite  expan¬ 
sions  w  hich  are  valid  throughout  the  region  of  interest. 

For  the  inner  expansion,  we  as-'ume  that  the  physical  parameters  R“  =  0  and  A/  =  e  are 
related  (dependent)  by  a  xet  of  mathematical  parameters:  gi.  C;-  4.1  •  ■  M-  by  a  "stretching 
transformation  of  the  form: 

£  =  eo + ,  (1-5) 

where  4  is  a  new  small  parameter  that  is  defined  by: 

5(|i)  =  (sgn5)4^  .  (1.6) 

The  positive  constants  a  and  b  are  to  be  determined  by  the  nature  of  the  singularities  of  F(X. 

£.  0)  near  e  =  eo,  where  Eq  is  the  singular  point  of  interest.  Let  the  dependent  variable 

z(4)  =  a(£(4),5(4))  be  written  as  the  expansion 

z  =  izj4J.  (1-7) 

1=0 

.Vote  that  for  the  pendulum  problem  F(a,  £.  0)  =  0  has  a  singular  point  at  £  =  £„  =  0.  The 
expansions  (1.7)  are  called  the  "inner  expansions"’  and  the  ZjS  are  called  the  "inner  coefficients 
Substituting  (1.5),  ( 1.6)  and  (1.7)  into  ( 1.2).  and  by  simplifying  the  inner  expansions  with  a=b=2 
one  obtains  the  following  solutions: 
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Zi  =  1  + 

Z2  =  I  + 


+  4 


I  - 


V^^  +  4  J 


-q-  Vq"  +4 


•  2 

.  .  2  ,  "5 

f 

+ 

\ _ 

^  T 

VV+4, 

^^+0(^1^) , 
|i^  +0(4^) . 


A/  =  e  =  qia^ 


(1.9) 


i.lO) 


and 

R-=5  =  |a-.  (1.11) 

Keeping  )i  fixed  and  taking  the  limit  |q|  ^  it  is  easy  to  see  that  zi  matches  asymptoti¬ 
cally  with  Xi  for  q  -» ■»  and  with  a.2  for  q  ^  -  =«.  Similarly  the  inner  solution  Z2  matches 
asymptotically  (|^|  — >  <»)  with  and  Xi-  New  combining  the  inner  and  the  outer  expansions 
appropriately,  we  get  the  composite  expansions: 


X,  = 


1  +  r2  +  (1  +  _L)R4 
A/ 


Al 


(l-u(A/))  + 


1  + 


2R2  -  a/  -  rW(A//R^)^  +4 


^  1  +  — ^ _ 

2  [  rW(A//r2)2-^4 

_L^  +  r2-(1  +  ±)r4 
1  A/  Al 


u(A/)-(l  -u(A/)) 


1  +  R-  +  (1  +  ^IR*^ 


-  u(A/) 


1  -Af  +A/2  +r2  -(1  +  — IR-^ 

Al 


+  0(R^) 


(i.i; 
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Xj  — 


(1  -u(Ai))  + 


1  +  2R^-A/+rW(A//R^)^  -^4 


A/Mi_ 


A/ 


rW(A//R^)-  +4 


+ 


1  +R^  +  (1  +  — jR-^  I  U(A/)-  !  l-A/+A/‘+R‘-(l  +  -^-jR-^  |{1  -  u(A/)) 
A/  I  I  A/ 


I 


-  [  1  +  R^  +  n  +  ^  iR^  |u( A/ )  +  0(R^  )  . 


where 


u(A/)  = 


1  .  A/>0 


0  ,  A/  <  0 

The  plots  of  eigenfrequencies  ,  Xi  versus  A/  are  given  in  Figure  3  for  both  the  exact  solu¬ 
tions  and  the  solutions  obtained  above  by  the  singular  perturbation  technique.  These  are  in 
excellent  agreement.  Thus,  the  singular  penurbation  technique  leads  to  qualitatively  correct 
asymptotic  approximations  that  are  often  close  to  true  solutions  and  can  be  used  as  a  mathemati¬ 
cal  tool  to  generate  quantitatively  accurate  solutions  for  a  wide  variety  of  linear  and  nonlinear 
structural  dynamics  problems.  The  methodology  is  general  and  systematic  and  when  combined 
with  elementary  singularity  theory,  should  provide  a  powerful  technique  to  study  the  mode  local¬ 
ization  phenomenon  in  any  finite  order  linear  or  cyclic  dynamic  chain. 
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ABSTRACT 

An  investigation  of  the  eigenvalue  loci  veering  and  mode  localization  phenomenon  is 
presented  for  mistuned  structural  systems.  Examples  from  both,  the  weakly  coupled  uniaxial 
component  systems  and  the  cyclic  symmetric  systems,  are  considered.  The  analysis  is  based  on 
the  singular  pemu’bation  techniques.  It  is  shown  that  uniform  asymptotic  expansions  for  the 
eigenvalues  and  eigenvectors  can  be  constructed  in  terms  of  the  mistuning  parameters  and  these 
solutions  are  in  excellent  agreement  with  the  exact  solutions.  The  asymptotic  expansions  are 
then  used  to  clearly  show  how  singular  behavior  in  the  eigenfunctions  or  modeshapes  leads  to 
mode  localization. 
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1.  INTRODUCTION 

In  recent  years  there  has  been  tremendous  interest  in  the  vibrations  and  structural  dynamics 
community  in  the  phenomenon  of  mode  localizadon.  This  interest  stems  fiom  the  recognition 
that  large  systems  composed  of  nominally  identical  subsystems  inevitably  involve  minor  devia¬ 
tions  from  the  idealized  structures  and  these  disorders  or  perturbations  can,  under  appropriate 
conditions,  cause  unexpectedly  large  deviations  from  the  predicted  behavior  in  the  nominal  or 
idealized  system  modes.  Important  technical  applications  where  these  problems  arise  include 
mistuned  bladed  disk  assemblies  [1,2]  and  large  space  structures  [3,4]. 

It  is  well  understood  by  now  that  the  presence  of  small  irregularities  in  nearly  periodic  struc¬ 
tures  may  inhibit  the  propagation  of  vibration  and  localize  the  vibration  modes.  Depending  on 
the  magnitude  of  perturbations  (disorder)  in  the  individual  components  and  the  strength  of  inter¬ 
nal  coupling  between  the  subsystems,  the  mode  shapes  may  undergo  dramatic  changes  and 
become  strongly  localized  when  small  perturbations  are  introduced,  thereby  confining  the  energy 
associated  with  a  given  mode  to  a  small  geometric  region.  This  phenomenon  is  referred  to  as 
mode  localization.  Pierre  [S]  suggested  that  strong  mode  localization  and  eigenvalue  curve  veer¬ 
ing  are  two  manifestations  of  the  same  phenomenon.  Therefore,  the  investigation  of  the  curves 
of  the  eigenvalues  or  natural  frequencies  in  the  neighborhood  of  the  ordered  state  is  sufficient  for 
determining  the  occurrence  of  strong  mode  localization.  That  eigenloci  veering  phenomenon 
can  occur  in  disordered  structures  under  certain  conditions  has  been  explained  qualitatively 
using  geometric  arguments  in  [6];  where  it  was  also  hinted  that  quantitative  results  might  be 
obtainable  by  the  use  of  singular  perturbation  analysis. 

For  mistuned  linear  chains,  Pierre  [5],  and  Pierre  and  Dowell  [7]  showed  that  the  straightfor- 
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ward  expansion  in  tenns  of  mistiming  parameters  breaksdown  in  the  case  of  weak  coupling. 
This  arises  because  the  idealized  system  that  is  being  perturbed  has  natural  frequencies  with 
multiplicity  >  1.  They  then  developed  a  so-called  "modified"  perturbation  technique  which  pro¬ 
vided  a  good  approximation  to  the  exact  eigenf^uencies  and  showed  good  agreement  with 
experimental  results.  In  the  case  of  strong  coupling  between  identical  subsystems,  no  such 
difficulty  arises  and  regular  penurbation  expansions  in  terms  of  mistuning  parameters  are  uni¬ 
formly  valid. 

For  systems  with  cyclic  symmetry  or  spatial  periodicity,  however,  mode  localization  can 
arise  in  the  presence  of  perturbations  which  split  the  degenerate  or  coincident  eigenvalues,^ 
irrespective  of  the  strength  of  internal  coupling  [6].  Using  differential  topological  ideas  it  was 
shown  qualitatively  in  [6]  that  circularly  configured  systems  which  have  cyclic  symmetry  exhibit 
complicated  topological  behavior  even  for  strong  coupling  when  small  perturbations  are 
imposed.  Furthermore,  the  frequency  response  of  a  perturbed  cyclic  system  depends 
significantly  on  the  form  of  the  perturbation.  Such  cyclic  periodic  systems  are  important  to  the 
analysis  of  vibrations  of  bladed  disk  assemblies. 

In  the  present  work  we  show  that  the  singularity  causing  the  breakdown  of  the  straightfor¬ 
ward  expansion  can  be  analyzed  by  the  well  developed  singular  perturbation  techniques  [8]  and 
appropriate  uniform  asymptotic  expansions  for  the  eigenfrequencies  and  eigenvectors  can  be 
constructed  which  provide  a  correct  qualitative  and  good  quantitative  approximation.  Prelim¬ 
inary  results  on  eigenvalue  veering  for  the  now  standard  example  [5]  of  the  coupled  penduli 
shown  in  Figure  1  were  recently  reported  in  a  short  paper  [9].  Here  we  present  complete  details 
of  the  singular  penurbation  analysis  for  the  eigenvalue  problem  of  the  coupled  penduli  system. 
Using  the  uniform  asymptotic  expansions  for  eigenvectors,  the  occurance  of  mode  localization  is 
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then  related  to  sensitivity  with  respect  to  parameter  variations.  We  then  consider  the  simplest  of 
examples  of  systems  with  cyclic  symmetry  consisting  of  three  identical  masses  arranged  in  a 
ring,  interconnected  by  identical  springs  and  having  individual  torsional  stiffnesses.  The  eigen¬ 
value  veering  is  here  shown  to  exist  even  for  the  strong  coupling  case.  Finally,  based  on  the 
solutions  for  the  strong  coupling  case,  behavior  for  the  weak  coupling  limit  is  explored. 

The  basic  idea  of  analysis  by  the  singular  perturbation  technique  is  the  following:  by  apply¬ 
ing  the  regular  perturbation  technique  to  the  eigenvalue  problem  A(|>  =  A.0,  of  the  system,  we  can 
obtain  algebraic  expressions  for  the  eigenvalues  and  eigenfuctions  as  a  power  series  in  the  small 
parameter  or  perturbation  (say  5).  The  coefficients  of  the  power  series  are  dependent  on  s~ 
second  parameter  e  and  these  expansions  are  valid  for  sufficiently  small  8,  for  all  values  of  e,  so 
long  as  no  singularities  arise.  Singularities  occur  for  values  of  e  where  the  eigenfrequencies  and 
eigenfunctions  lose  their  smoothness  and  it  is  said  that  the  expansion  is  not  uniformly  valid  for 
all  £.  Away  from  the  singular  values  of  the  parameter  (e),  the  straightforward  expansions  are 
good  approximations  and  are  called  the  "outer  expansions".  The  neighborhood  of  the  singular 
parameter  point  is  then  stretched  or  rescaled  in  terms  of  a  new  parameter  so  as  to  remove  the 
singularity.  The  expansion  in  terms  of  the  new  parameter  is  valid  only  in  the  neighborhood  of 
the  singular  point  and  is  called  the  "inner  expansion".  The  inner  and  the  outer  solutions  can  be 
matched  where  their  domains  of  validity  overlap  and  then  a  composite  expansion  can  be  con¬ 
structed  which  is  valid  uniformly  throughout  the  function  domain  for  all  values  of  the  parameter 
e. 
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2.  THE  COUPLED  PENDULI 
2.1  Singular  Perturbation  Analysis 

Consider  the  system  consisting  of  two  weakly  coupled  penduli  as  shown  in  Figure  1.  The 
two  imponant  parameters  are  the  dimensionless  length  Af,  and  the  dimensionless  coupling 
between  the  two  pendulums  5  =  =  (k/m)/(g//).  The  dimensionless  parameter  A/  represents  the 

disorder  or  perturbation  in  the  individual  pendulums.  It  is  important  to  point  out  that  the  dimen¬ 
sionless  coupling  between  the  two  pendulums  5  «  1  for  weak  coupling  irrespective  of  A/.  The 
resulting  eigenvalue  problem  in  symmetric  form  generated  by  the  above  system  is  given  by 

A0  =  X<j>,  (1) 

where  5  =  =  co^/coj  ,  ©i  =  k/m ,  ©J  =g// , 

\+R^  -r2 

A  = 

[-R2  ‘ 

For  small  values  of  5,  it  is  natural  to  expand  eigenvalues  and  eigenfuctions  in  the  regular  expan¬ 
sion  as  powers  of  5  regarding  A/  as  a  parameter  in  the  range  of  interest.  Thus  we  write  A,  X,  and 
4)  in  powers  of  5  as 


A  =  A<j  +  Ai5  A25^  +  0(8^) , 

(2) 

X  =  Xq  -t-  Xj  8  ^28^  -t-  0(8^) , 

(3) 

<j>  =  4>o  +  8  +  <1>25^  +  0(8^ ) . 

(4) 

Substituting  (2),  (3)  and  (4)  into  (1),  equating  coefficients  of  each  power  of  8  to  zero,  and  solv¬ 
ing  the  resulting  sequence  of  homogeneous  and  nonhomogeneous  linear  systems  gives  the  fol¬ 
lowing  expansions  for  the  eigenvalues  and  the  corresponding  eigenvectors 
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=  1  +5+  [(1  +A/)/A/j  52+0(5^) . 


X2  =  l/(l  +  A/)  +  5-  (l+A/)/A/j52+0(5^). 


<>1  = 


0 

(1  +A/)/A/ 


5  + 


52  +  0(5^), 


(5) 

(6) 

(7) 


02  = 


0 

Cl 


(1  +A/)/Ai 
0 


5  + 


52+0(5^) 


(8) 


In  (7)-(8),  C  and  Ci  are  arbitrary  constants.  The  expressions  (5)  -  (8)  are  the  regular  (outer) 
expansions  of  the  eigenvalue  problem  for  small  coupling  8  which  depend  on  the  parameter  A/. 
These  are  valid  for  sufficiendy  small  8  for  all  values  of  Al.  When  Ai  0,  (5)  -  (8)  become 
unbounded  and  the  continuity  of  eigenvalues  and  eigenvectors  with  respect  to  the  perturbation  Al 
breaks  down.  Thus,  in  the  neighborhood  of  Ai  =  0  ihe  expansions  (5)-(8)  become  nonuniform 
and  singular  or  non-analytic  points  have  therefore  been  identified.  These  eigenvalues  in  (5)  and 
(6)  are  plotted  in  Figure  2  for  some  small  but  fixed  6  as  a  function  of  Ai.  Each  eigencurve  has 
two  branches,  one  valid  for  A/  >  0  and  the  other  for  Al  <0.  These  branches  are  indicated  by  the 
superscripts  *+’  and  and  correspond  to  A/  >  0  and  Al  <0,  respectively.  Note  that  since  Al  and 
5  have  been  treated  as  two  independent  parameters  there  is  no  control  over  expressions  (5)  -  (8) 
in  the  limiting  process  when  A/  0  and  5  ->  0.  By  stretching  the  neighborhood  of  the  singular 
parameter  value  A/  =  0  and  by  taking  into  consideration  the  nature  of  the  singularity  we  can  find 
an  exact  relation  between  these  two  parameters.  Then  Al  and  5  become  dependent  in  the  neigh¬ 
borhood  of  the  singular  parameter  value,  called  the  "inner  region".  The  solutions  of  the  problem 
in  the  inner  region  are  called  the  inner  expansions.  Asymptotically  matching  the  inner  and  the 
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outer  expansions,  and  combining  them  appropriately  then  gives  the  composite  expansions  which 
are  valid  throughout  the  interval  of  interest  in  A/  for  sufficiently  small  5.  Before  finding  inner 
expansions  for  the  eigenvalues  and  the  eigenfunctions,  we  mass  normalize  eigenfunctions  0^  and 
0^  to  get  a  unique  set  of  eigenfunctions 


In  the  inner  expansion,  we  assume  that  the  physical  parameters  =  8  and  A/  =  e  are  related 
(dependent)  by  a  set  of  mathematical  parameters  ^2»  ^3»  •  •  •  and  p  through  a  ''stretching" 
transformation  of  the  form 

e  =  eo+^li»  +  i^j(li‘y .  (11) 

>=o 

where  p  is  a  new  small  parameter  that  is  defined  by: 

5(^)  =  (sgn5)^^  (12) 

For  fixed  p,  the  quantity  ^  serves  as  the  internal  variable.  The  positive  constants  a  and  b  are  to 
be  determined  by  the  nature  of  the  singularities  of  the  characteristic  equation  F(X,e,5)  =  0  of  (1) 
near  e  =  Eq,  where  Eq  is  the  singular  point  of  interest.  Let  the  eigenfunctions  0  and  the  eigen¬ 
values  X  be  written,  in  the  inner  region,  as  the  expansions 
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0(e,5)  =  z(e(n),6(4))  *  z(n)  =  £  Zjpj ,  (13) 

j=o 

X(e,  5)  =  X(e(ii).  5(h))  =  Q  (h)  =  £  .  (14) 

>=o 

Note  that  for  the  pendulum  problem  F(X,e.O)  =  0  has  a  singular  point  at  e  =  e©  =  0.  The  expan¬ 
sions  (13)  and  (14)  are  the  "inner  expansions"  and  the  Zj’s  and  Qj’s  are  called  the  "inner 
coefficients".  Substituting  (11),  (12),  (13),  and  (14)  into  (1)  simplifying  the  inner  expansions 
with  a  =  b  =  I,  and  solving  the  sequence  of  eigenvalue  problems  obtained  by  equating  each 
power  of  H  to  zero,  we  obtain  the  following  inner  eigenvalues  and  inner  eigenfunctions 

z'  and  z^ 


(19) 


-lo¬ 


rn  = 


2 


(20) 


ao  = 


1 


V 


l+k^ 


12 


+  (1+^11)^ 


and 


b„  = 


V 


m  + 


1  + 


12 


+  (1  +  ^11)^ 


(21) 


(22) 


Here  the  eigenvectors  and  have  been  mass  normalized.  Keeping  |i  fixed  and  taking  the 
limit  1  ^  ->  <»,  it  is  easy  to  see  that  matches  asymptotically  with  Xj  for  ^  and  with  Xl 
for  ^  -«>.  Similarly  the  inner  solution  matches  asymptotically  ( (^  -> «»)  with  Xl  and  X|. 

In  fact,  it  can  be  easily  shown  that  the  outer  eigenvectore  and  (j>m  match  with  the  inner  eigen¬ 
vectors  z^  and  z^  in  exactly  the  way  the  eigenvalue  branches  match.  The  composite  expansions 
for  the  eigenvalues  and  the  eigenfunctions  are  now  obtained  by  combining  the  inner  and  the 
outer  expansions,  and  subtracting  the  common  part  of  the  two  expansions.  The  resulting  expres¬ 
sions  for  the  eigenvalues  and  the  eigenvectors  are 


X 


1  _ 
comp  “ 


1 

1-t-A/ 


-»-5- 


" 

■ 

r  ^ 

AI 

52 

(l-u(Ai))-^ 

l-t-5-t- 

u-L 

A/ 

52 

a;  j 

L  j 

u(A/) 


^  % 

• 

(A//6)2 

1-1- 

2-(A//5)-f-V(A//5)2-»-4 

2 

5-h 

\  + 

(A//5)-»-V(A//5)2-h4 

2 

^  4 

T 

6^ 
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-(1-u(A/))(1-A/+A/2+8-(1+1/A/)5^)-(1  +  6+(1  +  1/A/)62)u(A/)+0(6^),  (23) 


= 

^comp 


1  +  5+ 


1  + 


A/ 


(l-u(Ai))  + 


1 


l+A/ 


+  5- 


5* 


u(A/) 


r  1 - ; — 1 

r 

(A//5)2 

+  * 

1  + 

2-(A//5)-V(A//5)2+4 

8+ 

1  + 

*  ^ 
(A//5)  +  V(Ai/5)2+4 

T 

52 

L  J 

L  2  J 

J 

-(1-u(A/))(1+5+(1  +  -^)82)-u(A/)[1 -A/ +  A/2  +  6-(1 +  47)5^1  +  0(5^)*  (2-^) 

A/  A/ 


-1+5/2A/ 

<t)  W  +  u(A/)<|>iL  +  (l-u(A/))({>^  -  u(A/) 

8/A/--^ 

+ 

A/ 

.  582 
”2Ai 

7AP 

l-u(A/) 


A/ 


1  + 


82 


2A/^ 


1  + 


82 


2A/- 


0 

-A/ 


8(1  +  -^) 

2A/2 


0(82) , 


which  simplifies  to 


<t> 


comp 


=  ao1 


p  1 

k2A/2/62 

-k 

1 

+ 

l+k2 

8 

1 

0 

d 

8  U  0(82), 


(25) 
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and 

•  « 

r_,i  fm^A/Wl 

<t>comp  =  boj  J  +  1+m^  5  ^+0(5^),  (26) 

M-  -1  0 


where  k,  m,  a^,,  and  bo  are  already  defined  in  equations  (19),  (20),  (21)  and  (22)  respectively  and, 


u(Al)  = 


1  ,  A/  >  0 
0,  A/<0‘ 


The  composite  solutions  (23)  -  (26)  of  the  eigenvalue  problem  (1)  are  the  asymptotic  approxima-  - 
tions  for  small  coupling  5,  and  are  uniformly  valid  for  all  mistunings  A/. 

The  eigenfiequencies  and  X,|)mp,  obtained  by  the  singular  perturbation  analysis,  arc 
plotted  in  Figure  3  as  a  function  of  the  parameter  A/.  Since  the  coupled  penduli  system  is  sim¬ 
ple,  exact  expressions  for  the  eigenfiequencies  are  easily  obtained  and  they  are  also  plotted  in 
the  figure.  The  exact  solutions  for  the  eigenvalue  problem  are  given  in  the  Appendix.  Clearly, 
there  is  an  excellent  agreement  in  the  exact  frequencies  and  their  asymptotic  approximations. 
The  asymptotic  solutions  also  clearly  display  the  veering  phenomenon. 

We  now  study  the  behavior  of  eigenvectors  for  the  case  of  weak  coupling  when  A/  -  0. 


2.2  Eigenvector  Rotations  and  the  Sensitivity  Function 

In  our  earlier  work  [6]  with  mistuned  cyclic  systems,  it  was  suggested  that  localization  of 
modes  can  be  investigated  by  considering  the  sensitivity  function,  and  the  rotations  of  eigenvec¬ 
tors  under  variations  of  parameters.  The  sensitivity  function  of  eigenvectors  or  eigenvector  sen- 
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sitivity,  in  short,  is  defined  as 

1 1  Su  1 1  =  ^/^(SuSu)  (Frobenius  norm) , 

where  Sy  =  Uq^Au,  Su  denotes  the  complex  conjugate  of  transpose  of  Su,  Uq  is  the  modal  matrix 
of  eigenvectors  for  zero  coupling  (5  =  0),  ‘tr’  denotes  the  trace  of  the  matrix,  and  Au  =  u  -  u^ 
where  u  is  the  modal  matrix  for  nonzero  coupling  (S  ^  0).  The  eigenvector  sensitivity  evaluated 
for  the  pendulum  problem  turns  out  to  be 

I  I  Su  I  I  =  — p=-  "n/C-Pii  +P2i)^  +  (Pi1  +F2i)^  +  (“Pi2+P22)^  +  (Pi2+P22)^  .  (27) 

V2 

where 


^  ,  1  aok2(A/)2 

Pll  =-aok+  — zsr  + 


^J2  5(1 +k^) 


1 

Pi2  =  bom - =•  + 


^j2  5(1  +m^)  ’ 


P21  ~  ^  I —  » 

P22  =  bo - ]=■ . 

V2 

We  can  also  define  the  angles  between  the  eigenvectors  and  <l)o  for  the  unperturbed  (5  =  0) 
system,  and  the  eigenvectors  <l>comp  for  tbe  penurbed  system.  These  angles  are  given 

by 


«t>o  _  ^-Pl 

II  00 II II  0eomp  II  VT'^Jl  +  Pf 


(28) 
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^  »  ^oamp  ^  ^  P2 

old'll  0l.mpll  " 

where 


Pi  =k  + 


k^(A/)^ 
dd+k^)  ’ 


P2  =  m  + 


m^(A/)^ 

6(l+ni^) 


Plots  of  the  sensitivity  function  of  eigenvectors  and  the  cosine  of  the  angle  between  nominal 
(5  =  0)  and  perturbed  eigenvectors,  as  a  function  of  the  mistuning  A/,  are  given  in  Figures  4  and 
5,  for  both  the  exact  solutions  and  the  asymptotic  approximations  obtained  above  by  the  singular 
perturbation  technique.  These  solutions  are  in  excellent  agreement  Figures  4  and  5  clearly 
show  and  confirm  the  expectation  that  the  eigenvectors  for  the  weakly  coupled  system  undergo 
rapid  changes  in  the  vicinity  of  the  singular  point  Furthermore,  either  of  the  two  criterion  can 
be  effectively  used  as  a  quantitative  measure  and  indicator  of  the  mode  localization 
phenomenon. 


3  CYCLIC  SYSTEMS 

Consider  three  identical  particles,  each  of  mass  m  arranged  in  a  ring  and  interconnected  by 
identical  springs  of  stiffness  k^.  Assume  that  all  the  masses  are  hinged  to  the  ground  by  torsional 
springs  of  stiffness  k,  and  that  the  radius  of  the  ring  is  r,  as  shown  in  Figure  6.  As  a  perturbed 
system  we  consider  the  case  when  two  of  the  torsional  springs  are  perturbed  by  £1  and  £2-  The 
eigenvalue  problem  corresponding  to  this  system  is  given  by 
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a+Ei  -<o|  -<oi 

-0)|  a+e2  -cal  0  =  (29) 

-co|  -<i)^  a 

2kc  kf  2  kc 

where  a  = - h - ^  and  (Oc  =  —  • 

m  mr^  m 

As  was  the  case  for  the  example  of  linear  chains,  the  interest  here  is  in  the  development  of 
asymptotic  expansions  for  the  eigenvalues  and  eigenvectors  in  terms  of  the  perturbation  parame¬ 
ters  El  and  £2.  As  we  shall  see  the  unperturbed  cyclic  system  (£1  =  £2  =  0),  has  a  double  eigen¬ 
value  and  an  isolated  eigenvalue,  and  thus,  introduction  of  perturbations  is  expected  to  split  the 
degenerate  eigenvalue  pair.  Expecting  that  the  eigenvector  behavior  will  be  governed  by  the 
eigenvalue  behavior,  similar  to  the  case  of  the  coupled  penduli,  we  restrict  the  developments  to 
expansions  for  only  the  eigenvalues  of  the  perturbed  three  particle  system. 

The  eigenvalue  problem  (29)  results  in  the  characteristic  equation 

F(X,,£i  ,£2)  =  — (3a  +  £i  -♦■£2)^^  +  (3a^  -♦■2a£2  +  2a£2  +61C2  ~  30)^)^ 

4-  2o)c  +  cOcOa  £1  +  £2)  -  (a^  +  a£i)(a  +  £2) j  =  0 .  (30) 

First  consider  the  unperturbed  cyclic  system  the  eigenvalues  or  roots  of  (30)  are  given  by 
=  a  -4  cOc  .  =  a  -  2cDc  • 

The  corresponding  roots  of  (31)  are 

—1—2  2  ->2 
X  =X  ,  X  =-2to‘ . 

Thus,  there  is  a  coincident  pair  of  eigenvalues  and  one  isolated  eigenvalue.  So  long  as 
(Oc  “  0(1),  the  two  distinct  eigenvalue  arc  well  separated. 
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To  study  the  perturbed  problem  we  first  introduce  the  coordinate  transformation  x  =  X.  —  a,  so 
that  equation  (30)  results  in 

-  (ei+€2)X^  +  (eie2  “  3t0c)X  +  o)^(2a)|  +  £+£2)  =  0  (31) 

We  can  express  the  solutions  to  (31)  as  regular  functions  of  the  parameters  £},  £2,  and  cOg  by 

00 

writing  x(£i.£.)  =  2Xj(£i)£i  as  a  power  series  in  £2.  Substituting  the  resulting  expression  in 
i=o 

(31)  and  proceeding  in  the  usual  manner,  the  expansions  for  the  three  roots  of  (31)  turn  out  to  be 


Wc-Ei 


a)|£i 


£2+0(62). 


(32) 


t 


X^  =  X02  +  Xl2e2  +  Xl2^\  +  0(£^)  , 


X^=Xo3+Xi3e2+X23e^  +  0(£i), 


where 


_  £1  -  0)2  +  Pi  Xra-Xo2ei  -  Wc 

X02  ^  .  X12  _  2  /»  _ 

^  3xo2~2xo2£i  “  3o); 


X22  = 


Xi2(ei-3xo2)  +  Xi2(2Xo2-£i) 

3x^2-2X02ei  -  3o)^ 


_  £1 -0)2  -  Pi  Xra-XoaEi  - 

X03 - -I -  .  Xl3  = 


3x^3-2X0361  -  3o)^  ’ 


X23  = 


X  13(61-3x02) +  Xi3(2X03-6i) 
3x^3-2X0361  -  30)^ 


and 


(33) 

(34) 
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Pi  =  +  2c(>?ei  +  e?  . 

The  expressions  (32)  -  (34)  are  the  straightforward  expansions  of  the  eigenvalue  problem  for 
small  £2-  When  ei  ->0,  yi}  and  become  unbounded  and  the  continuity  of  the  eigenvalues 
with  respect  to  the  parameter  or  perturbation  Ej  breaks  down.  The  third  eigenvalue  x^  always 
remain  bounded  and  continuity  is  preserved  for  all  values  of  Ei  so  long  as  the  interconnecting  or 
coupling  spring  constant  kc  is  0(1).  These  eigenvalues  (32)-(34)  are  shown  in  Figure  7.  Clearly 
El  =  0,  E2  =  0  is  a  singular  point  of  the  expansions  and  (32)-(34)  are  the  outer  expansions  valid 
for  small  E2  away  from  £1  =  0. 

Inner  expansions,  which  are  valid  in  the  neighborhood  of  singular  point  or  the  parameter 
values  where  the  outer  expansions  breakdown  are  now  obtained.  The  expansion  process  for  the 
cyclic  system  is  very  similar  to  the  one  presented  in  section  2  for  the  linear  chain  system.  Thus, 
the  physical  parameter  perturbations  £1  and  £2  are  related  to  a  set  of  parameters  ,  ^2>  ^3>  ti 
via  a  "stretching"  transformation  of  the  form 

£  =  +  (35) 

j=o 

where  p.  is  a  new  small  parameters  that  is  defined  by 

6(p)  =  (sgn  6)n*’ .  (36) 

The  positive  constants  a  and  b  are  to  be  determined  by  the  nature  of  the  singularities  of 

F(X,£i,£2)  =  0  near  Ej  =0,  the  singular  point  of  interest.  Let  the  dependent  variable  or  the 

eigenvalue  Q(|i)  =  fl(£i(n),  £2(p))  be  written  as  the  expansion 

n  =  £njnj. 

j=o 


(37) 
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Note  that  for  the  cyclic  system  F(X,ei.e2  =0)  has  a  singular  point  at  ei  =  £2  =  0.  The  expan¬ 
sions  (37)  are  called  the  "inner  expressions"  and  the  QjS  are  called  the  "inner  coefficients".  Sub¬ 
stituting  (35),  (36)  and  (37)  into  (30),  simplifying  the  inner  expansions  with  a=b=l,  and  perform¬ 
ing  the  permrbation  analysis,  the  following  roots  of  (30)  are  obtained 


=  tOc  +  flu  |1  4- 


(4+i)flii  -flii-^flii 


2a)iV4M+l 


fl^  =  col  fli2  ^l-J- 


fl^  =-20)2-4 
where 


4+1 


11  + 


V*  J 


-2 


270)2 


(^2-^-H)^2  +  0(^3)  , 


(38) 


(39) 


(40) 


El  =  ,  62  =  H  ,  All  =  Q  ^  ^-fl) 

3  3 

Keeping  fixed  and  taking  the  limit  |  ^  —>  00,  it  is  easy  to  see  that  fl^  matches  asymptotically 
with  for  ^  — 4  00  and  with  xi  for  ^  — »  —00.  Similarly  the  inner  solution  fl^  matches  asymptoti¬ 
cally  (1^  ->  00)  with  xi  and  X-.  As  expected,  the  third  root  fl^  automatically  matches  with  x^ 
as  no  singular  behavior  is  displayed  in  this  case.  Now  combining  the  inner  and  the  outer  expan¬ 
sions  appropriately,  we  get  the  composite  expansions 

Xcomp  =  X^(u(Ci ))  +  x'  ( l“U(ei ))  +  fl*  -  u(ei )  (common  parts  of  x^  ,  fl* ) 


-  (l-u(ei))  (common  parts  of  x*  .  A*) , 
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xlomp  =  X*(u(ei))  +  X^(l-u(ei))  +  -  u(ei)  (common  parts  of 

-  (l-u(ei))(common  parts  of  . 


Xcomp  —  X  > 


where  u(ei )  = 


1  ,  ei  ^0 

|0  ,  ei  <0  • 


Transforming  these  expansions  back  into  the  original  coordinates  gives  the  following  composite 
expansions  for  eigenvalues 

^comp  ~  ^11  ^2 


(1  +  ei/e2)i2ii  -  flu  -  (£1/62)^11 


^comp  ~  +  i2i2  £2 


Qi2  +  £1/82^12  ~  (£i/£2  1)^12 


ei  +  0(ei) , 


(41) 


2coi^(ei/e2)^-(ei/e2)+l 


ei  +  0(ei) , 


^»mp  =  X03  +  (Xl3)e2  +  X23(ei)  +  0(£i)  , 
where 


(42) 


(43) 


«ii  = 


(£1/62  +  1)  +  V(^i/£2)^  -(ei/£2)+  I 


ni2  - 


(ei/£2  +  1)  -  y/(ei/£2)^-(ei/e2)+ I 
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-XoBgi 

3xfe  -  2X0361  -  3(0^  ’ 

El  -  (0^  -  Pi 

Xo3  = - 2  ’ 

Pi  =  >/9tOc  +  2co^ei  +  El  . 

The  plots  of  eigenfirequencies  Xcomp.  ^comp  and  X^omp  as  a  function  of  Ei  are  given  in  Figure  8a 
and  8  b.  The  third  eigenvalue  X^^p  always  behaves  as  a  regular  function  as  is  clear  from  Figure 
8a.  These  results  show  that  for  the  cyclic  system  with  strong  coupling  the  effect  of  perturbation 
is  identical  to  that  for  the  pendulum  problem  with  weak  coupling.  That  is,  the  curve  veering  of 
the  eigenvalues  X^omp  and  X^omp  is  almost  the  same  as  the  curve  veering  of  the  weakly  coupled 
pendulum.  Consequently,  it  is  expected  that  the  mode  localization  of  the  eigenvectors  within  the 
strongly  coupled  cyclic  system  be  similar  to  that  in  the  weakly  coupled  pendulum  problem. 
Therefore,  the  results  should  be  obtained  for  the  singular  behavior  of  the  eigenvectors  of  the 
weakly  coupled  pendulum  problem  in  section  2  should  be  valid  quatitatively  for  the  cyclic  sys¬ 
tem. 

The  coupling  constant  k^  plays  a  very  important  role  in  the  eigenvalue  veering  behavior.  For 
the  cyclic  system,  kc  =  0  leads  to  three  coincident  eigenvalues  as  opposed  to  the  case  of  strong 
coupling  when  only  a  double  eigenvalue  appears.  The  composite  expansions  obtained  earlier 
were  determined  under  the  assumption  that  Iq;  -  0(1).  We  now  use  these  expansions  and 
explore  the  behavior  of  eigenvalues  in  the  limiting  ease  of  kc  — >  0.  Note  that  the  coupling  con¬ 
stant  oig  appear  in  the  denominator  of  the  asymptotic  expansions  (41)-(43).  As  coc  is  reduced 
the  isolated  eigenvalue  Xco^p  moves  closer  to  X^^p  the  variation  of  X  s  with  E2 
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becomes  very  rapid,  as  is  evident  from  Figures  9  and  10.  Figure  1 1  shows  the  three  composite 
eigenvalue  expansions  for  small  Ej  and  £2  as  a  function  of  kc  (or  (ol)  and  the  expansions  clearly 
breakdown  in  the  vicinity  of  k<;  =G.  Thus,  the  composite  expansions  (4l)-(43)  are  behaving  as 
outer  expansions  in  the  weak  coupling  limit.  Ft  should  be  possible  to  now  construct  uniformly 
valid  expansions  for  eigenvalues  as  a  function  of  the  coupling  constant  k^  by  one  more  use  of  the 
singular  perturbation  technique  whereby  the  neighborhood  of  Iq.  =  0  is  stretched  and  an  inner 
expansion  is  developed.  The  localization  and  veering  behavior  of  thus  obtained  composite 
eigenvalue  expansions,  which  will  be  valid  for  small  enough  £2  for  all  £1  and  Iq.,  is  expected  to 
be  much  more  interesting  and  is  being  presently  studied. 

4.  SUMMARY  AND  CONCLUSIONS 

Singular  perturbation  technique  has  been  applied  to  two  parameter  eigenvalue  problems  to 
obtain  uniformly  valid  algebraic  expansions  for  the  eigenvalues  and  the  eigenvectors  for  two 
example  systems.  Utilizing  these  expansions,  eigenloci  veering  and  the  mode  localization 
phenomenon  have  been  studied.  A  sensitivity  function  and  the  rotation  of  eigenvectors  have 
been  introduced  as  criteria  to  visualize  mode  localization  phenomenon  in  the  vicinity  of  singular 
points.  One  example,  that  of  the  two  weakly  coupled  penduli,  represents  systems  of  the  linear 
chain  type  with  only  one  weak  coupling  spring.  The  example  of  three  mass  particles  belongs  to 
strongly  coupled  systems  with  cyclic  symmetry.  For  the  coupled  penduli  system,  the  eigenvalue 
and  eigenvector  expansions  arc  found  to  be  in  excellend  agreement  with  the  exact  results. 

It  is  shown  that  eigenvalue  curve  veering  occurs  both  in  the  weakly  coupled  penduli  and  the 
strongly  coupled  cyclic  system.  The  effects  of  mistuning  perturbations  which  split  a  pair  of 
coincident  eigenvalues  is  identical  in  both  the  cases.  T1  e  eigenvector  sensitivity  function  and 


the  angle  of  rotation  of  eigenvectors  are  shown  to  be  two  equally  good  candidates  for  visualizing 
mode  localization  phenomenon  near  singular  points.  The  composite  expansions  for  the  per¬ 
turbed  cyclic  system,  which  are  uniformly  valid  in  the  ca.se  of  strong  coupling,  are  shown  to 
breakdown  in  the  limiting  case  when  the  coupling  stiffness  goes  to  zero.  This  clearly  is  related 
to  the  fact  that  all  the  eigenvalues  for  the  cyclic  system  in  the  weak  coupling  limit  are  identical. 
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APPENDIX 

It  can  be  easily  shown  that  the  two  free  vibration  natural  frequencies  on}  and  the  correspond¬ 
ing  mass  normalized  eigenvectors  x‘  for  the  coupled  pendula  problem,  obtained  from  the  exact 
solution  of  the  eigenvalue  problem,  arc  given  by 


2a)?,2  =2R^  +  \  + 


1  - 
+ 


1-t-A/ 


4R^  -I- 1  -  -L-  + - L- 

1+A/  (1+A/)2 


‘/4 
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The  sensitivity  of  eigenvectors  is  then  given  by 


IISTlI  =  -|=-  +  (qil+<l2l)^  +  (-<112+022)^  +  (qi2+q22)^  . 

V2 


where 


1  1 

qii  =aXii  +^. 


qi2  =  Pxii  ~  * 


q.=a-^. 


q22  -  p  -  , 


The  angles  between  the  nominal  and  the  perturbed  eigenvectors  are  then 


cos 


eT  = 


i-xli 


>^Vi+0cii 


i-x?i 

>^Vl+(Xn)^ 
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Figure  I .  Two  coupled  oscillators. 


’Comp 


Figure  8a.  Eigenvdues  from  composite  expansions  for  the  perturbed  cyclic  system;  kc  =  2, 
kt  =  1,  £2  =0.1,  m  =  1,  r  =  1. 
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Figure  8b.  Behavior  of  the  two  eigenvalues  X^p,  showing  curve  veering  for  the  strong 
coupling  case;  k<.  =  2,  kt  =  1,  £2  =  0.1,  m  =  1,  r  =  1. 


duio 


■comp 


Figure  8b.  Behavior  of  the  two  eigenvalues  Xlamp^  ^wmp  showing  curve  veering  for  the  strong 
coupling  case;  k,.  =  2,  kt  =  1,  £2  =  0.1,  m  =  1,  r  =  1, 
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Figure  10.  Behavior  of  the  composite  eigenvalues  as  a  function  of  the  coupling  parameter  kc; 
kt  =  l,ei  =0.02,ei  =  0.1,m  =  l.r=  1. 
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I.  Introduction 

An  imponant  subject  in  the  dynamics  and  control  of  structural  systems  is  the 
behavior  of  structures  under  transient  or  steady  state  excitations.  In  this  work,  we 
examine  the  stability  of  the  geometric  form  of  the  spatial  configuration  of  struc¬ 
tural  systems  when  the  structural  parameters  are  subject  to  small  perturbations, 
and  the  implications  of  this  instability  for  frequency  response.  We  show  that  cir¬ 
cularly  configured  systems  which  nominally  have  cyclic  symmetry  exhibit  com¬ 
plicated  topological  behavior  when  small  perturbations  are  impressed  on  them. 
We  further  show  that  the  frequency  response  of  a  perturbed  cyclic  system 

,  ‘  I’yKt  K..  ^  ?-•"  - 

'  -epr -c,..:  ,.'1  ,r.  • 
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depends  considerably  on  the  form  of  perturbation.  On  the  other  hand,  a  recdiinear 
configuration  of  nearly  identical  subsystems  does  not  ehibit  modal  insubility. 
Usually,  both  kinds  of  systems  are  implicidy  assumed  to  undergo  similar  qualita- 
five  behavior  under  a  small  perturbation  whereas,  in  fact,  the  cyclic  configuration 
exhibits  a  very  stange  behavior,  [1]. 

The  disnnciion  between  the  behavior  of  cyclic  and  rectilinear  configurations 
under  a  perturbation  is  important  because  many  engineering  structures  are  com¬ 
posed  of  identical  substructures  which  are  replicated  either  in  a  uni-axial  chain, 
or  in  a  closed  cyclic  formation  where  modal  control  is  of  interest.  Examples  of  the 
former  case  of  periodicity  occur  in  structures  such  as  space  platforms  and  bridges, 
which  have  an  obvious  periodicity  of  the  uni-axial  kind.  An  example  of  cyclic 
periodic  systems  is  a  turbine  rotor,  which  consists  of  a  set  of  nominally  identical 
blades  mounted  on  a  central  hub,  and  often  referred  to  as  a  “bladed  disk  assem¬ 
bly”  [2].  When  all  the  blades  are  truly  identical,  then  the  system  is  referred  to  in 
the  literature  as  a  tuned  bladed  disk  assembly.  Pracdcal  realities  of  manufacturing 
processes  preclude  the  existence  of  exact  uniformity  among  all  the  blades.  When 
residual  differences  from  one  blade  to  another — no  matter  how  small — are 
accounted  for  in  the  theoretical  model,  the  assembly  is  then  termed  a  mistuned 
bladed  disk. 

Our  primary  focus  in  this  investigation  is  on  bladed  disk  assemblies.  However, 
since  we  approach  the  problem  fiom  a  generalized  viewpoint,  the  conclusions  to 
be  drawn  will  be  of  relevance  to  other  periodic  systems.  Therefore,  in  the  sequel, 
we  borrow  the  ‘tuned’  and  ‘mistuned’  terminology  from  the  bladed  disk  litera¬ 
ture,  and  apply  it  to  repedove  systems  having  cyclic  or  uniaxial  periodicity.  Thus, 
in  a  tuned  periodic  system,  the  nominal  periodicity  is  preserved,  whereas  it  is  des¬ 
troyed  in  a  mistuned  system. 

If  we  examine  the  system  matrices  of  the  linear  and  cyclic  chains,  we  observe 
a  fundamental  difference  in  forms.  The  dynamical  matrix  of  the  linear  chain  is 
usually  banded.  Banded  matrices  are  fiequendy  encountered  in  structural  dynam¬ 
ics.  A  special  form  of  banded  matrices  that  is  of  interest  here  is  the  tri-diagonal 
form  aj)sO,|i-y  |>  1.  On  the  other  hand,  the  system  matrix  of  a  cyclic  chain  has  a 
circulant  submatrix,  or  is  entirely  cireulant  or  block  circulant  [3].  Circulant 
matrices  usually  arise  in  the  study  of  circular  systems.  They  have  interesting  pro¬ 
perties  that  set  them  apan  from  matrices  of  other  forms  [4].  We  note  that  all 
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ciiculants  commute  under  muldplicadon,  and  are  diagonalizable  by  the  fouiier 
matrix.  One  of  the  most  imponam  consequences  of  the  foregoing  is  that  the  cyclic 
chain  has  a  series  of  degenerate  eigenvalues,  whereas  the  eigenvalues  of  the  uni¬ 
axial  chain  are  all  simple. 

We  know  that  a  tuned  circuiant  matrix,  having  a  muldtude  of  degenerate 
eigenvalues,  lies  on  a  bifurcadon  set  [5].  Thus,  the  reducdon  of  such  matrices  to 
Jordan  normal  form  is  an  unstable  operadon  [6].  Consequendy,  if  a  non-singular 
deformadon  due  to  mistuning  is  applied  to  a  circuiant  matrix,  then  some  of  the 
eigenvectors  will  undergo  rapid  re-alignmenu  if  the  mistuning  leads  to  a  crossing 
of  the  bifurcation  set.  If  however,  no  crossing  of  the  bifurcadon  set  takes  place, 
then  the  tuned  system’s  eigenvectors  will  be  very  sublc,  preserving  their  align¬ 
ment  under  small  penurbadons.  In  contrast,  the  eigenvectors  of  a  tuned  banded 
matrix,  being  analytically  dependent  on  parameters,  are  not  generally  disoriented 
by  mistuning  until  the  eigenvalues  are  pathologically  close  [7]. 

If  one  examines  the  literature  in  structural  dynamics,  it  is  observed  that  some 
unusual  behavior  has  been  reponed  in  the  study  of  perturbed  cyclic  systems.  This 
has  been  the  case  in  various  studies  of  rings  [8],  circular  saws  [9],  and  other 
cyclic  structures  [10].  But  that  such  anomalous  behavior  is  due  to  a  “geometric 
instability"  inherent  in  the  cyclicity  of  the  tuned  system  has  not  been  previously 
established  in  the  literature,  to  our  knowledge.  Indeed,  it  is  often  assumed  (see, 
for  instance,  [11])  that  the  linear  and  cyclic  chains  would  undergo  the  same  quali- 
utive  behavior  under  slight  parameter  penurbadons  so  that  small  order  pemuta- 
rions  of  the  system  matrix  will  lead  to  no  more  than  small  order  differences  in  the 
system  response  relative  to  the  unperturbed  case,  if  the  system  has  “strong  cou¬ 
pling". 

In  this  paper,  we  show  that  such  an  assumption  regarding  qualiudve  behavior 
docs  not  actually  hold  in  the  case  of  cyclic  systems;  that  cyclic  systems  exhibit  a 
peculiarity  of  their  own  under  parameter  perturbation;  that,  although  a  certain 
amount  of  mistuning  may  produce  litde  diference  relative  to  the  tuned  datum  in 
one  case,  a  considerable  change  could  be  induced  if  a  slightly  different  kind  of 
mistuning  is  applied  to  the  same  cyclic  system;  that  such  apparently  erratic 
behavior  arises  in  cyclic  system,  even  when  the  system  has  “strong”  coupling.  In 
carrying  out  this  work,  we  borrow  from  certain  developments  in  differential  topol¬ 
ogy  specifically,  from  Arnold’s  monumental  work  in  singularity  theory  [6,  12-16]. 
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il.  Topological  Dynamica  of  Quadratic  Syatama 

In  mistuned  dynamical  systems,  a  major  concern  is  to  understand  which 
specific  kinds  of  mistuning  parameters,  or  combinations  thereof,  lead  to  unaccept¬ 
ably  high  amplitude  ratios.  In  this  section,  we  give  an  indication  of  the  taxonomy 
of  the  different  consequences  of  mistuning  in  the  hope  of  isolating  those  that  lead 
to  high  ratios. 

Consider  a  mechanical  system  under  small  oscillations  with  kinetic  and  poten¬ 
tial  energies  given  by: 

T  =  Ax'  Mx  >0,  U  =  'Ax'  Kx  >0-,  x,x*0.  (2.1) 

Under  the  influence  of  a  forcing  function  /  (r).  (2.1)  produces  the  following  equa¬ 
tions  of  motion  by  application  of  Lagrange’s  formula: 

Mx -^Kx  =/■  x.feCr  (2.2) 

where  M  and  K  are  symmetric  -'d  positive  definite.  A  theorem  in  linear 
algebra  shows  that  there  exists  some  non-singular  transformation  matrix  P  such 
that: 

P^'^P  =  /,  and  P^KP  =  A  (2.3) 

where  A  is  a  diagonal  matrix  of  eigenvalues  whose  elements  satisfy  the  equation: 

det(Af-Xi:)  =  0  (2.4) 

Consequently,  by  putting 

X  =  Pq.  (2.5) 

substituting  for  q  in  (2.1),  and  premultiplying  every  term  of  the  resultant  equation 
by  P^ ,  we  obtain  a  new  equation  set: 

q*  \q=  /',  (2.6) 

where  f'  =  P^f.  Hence: 

q,  +  X.q,  =  /',  for  /  =  1,  2,  •  ■  • ,  n.  (2.7) 

Systems  which  can  be  reduced  to  the  above  form  are  called  quadratic  systems. 
They  are  called  quadratic  cyclic  systems  if,  in  addition,  M  and  K  are  cyclic  or 
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cireulant  matrices.  Our  basic  aim  is  to  determine  the  nature  of  the  changes  in  the 
dynamical  properties  of  a  quadratic  system  of  a  given  order,  under  random 
diffimndal  perturbations  in  M  and/or  K.  Central  to  this  investigation  are  the  topo¬ 
logical  concepts  of  structural  stability  and  genericity. 

Let  A/  be  a  set  with  a  topology  and  an  equivalence  reladon  e.  An  element  x€N 
is  stable  (relative  to  e)  if  the  e-equivalence  class  of  x  contains  a  neighborhood  of 

X. 

A  property  P  of  elements  of  N  is  generic  if  the  set  of  all  Jt  6  N  satisfying  P 
contains  a  subset  A  which  is  a  counuble  intersection  of  open  dense  sets  [17]. 

Genericity  is  important  ^  our  context  because  a  generic  system  will  in  efBsct 
display  a  “typical"  behavior.  .More  concretely  if  a  given  generic  system  gives  a 
certain  frequency  response,  all  systems  produced  by  differential  parameter  pertur¬ 
bations  about  the  nominal  system  will  also  produce  frequency  response  curves 
that  are  not  only  slight  perturbanons  of  the  original  nominal  response  but  also 
geometrically  (isomorphic)  equivalent  to  it.  Such  systems  are  called  versal  defor¬ 
mations  of  the  nominal  system  [14],  A  versal  deformation  of  a  system  is  a  normal 
form  to  which  it  is  possible  to  reduce  not  only  a  suitable  representation  of  a  nomi¬ 
nal  system,  but  also  the  representation  of  all  nearby  systems  such  that  the  reduc- 
don  transformadon  depends  smoothly  on  parameters.  The  key  to  establishing  ver- 
sality,  and  hence  genericity,  is  the  topological  concept  of  transversality. 

Let  A/cM  be  a  smooth  submanifold  of  the  manifold  M.  Consider  a  smooth 
mapping  f  X  -» Af  of  the  parameter  space  P  into  A/;  and  let  d  be  a  point  in  F  such 
that  /  (11)  €  N. 

The  mapping  /  is  transversal  to  iV  at  n  if  the  tangent  space  to  .Vf  at  /  (ji)  is  the 
sum: 


7'Af/(n)  =  /• 

Consequently,  two  manifolds  intersect  transversally  if  either  they  do  not  intersect 
at  all  or  intersect  properly  such  that  perturbadons  of  the  manifolds  will  neither 
remove  the  intersecdon  nor  alter  the  type  of  intersecdon. 

Lemma  2.1,  see  ref  [14]. 

A  deformation  f(ii)  is  versal  if  and  only  if  the  mapping  fT  — »  M  is  transversal 
to  the  orbit  o/f  a;  ii  =  0. 
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The  abov*  result  is  crucially  important  because: 

(i)  It  classifies  from  the  set  of  all  perturbations  of  a  given  nominal  system, 
those  that  do  not  lead  to  radically  diftrent  dynamical  properties  from  the 
nominal. 

(ii)  It  separates  the  “good"  from  the  “bad”  perturbations  and  hence  enables  us 
to  concentrate  our  study  on  the  bad  penurbadons.  Let  Q  denote  the  family 
of  all  real  quadratic  systems  in  R".  The  set  Q  has  the  structure  of  a  vector 
space  of  dimension  Ainin  -r  1|).  It  can  be  shown  that  Q  also  has  the  struc- 
mre  of  a  differentiable  manifold  [13]. 

Let  Qv  denote  the  set  of  quadratic  systems  having  V2  eigenvalues  of 
multiplicity  2,  V3  eigenvalues  of  multiplicity  3  etc.  Qv  is  called  the  degen¬ 
erate  subfamily  of  Q. 


Theorem  2.1,  see  ref  [13]. 

The  transformation  fiF  ->  Q  is  transversal  to  Qv- 

Consequently,  a  generic  family  of  quadratic  sysi..ns  of  a  given  order  is  given 
by  a  transformadon,  /,  of  the  space  of  parameters  T  into  the  space  of  ail  quadratic 
systems  Q,  such  that  /  is  transversal  to  the  space  of  all  degenerate  quadranc  sys¬ 
tems  (2  v 

Hence  (2v  is  the  degenerate  (bad)  set  and  Q/Qv  ’s  the  generic  set.  Observe 
that  Q/Qv  and  are  transversal.  Consequently,  the  fundamental  group  of  the 
space  of  generic  real  quadranc  systems  is  isomorphic  to  the  manifold  of  systems 
without  degenerate  eigenvalues. 

The  above  discussion  leads  inevitably  to  the  following  conclusions; 

(i)  Radical  changes  in  the  dynamical  properdes  of  a  nominal  system  occurs 
under  perturbadons,  when  the  penurbadons  take  the  system  across  the 
boundary  from  QlQy  to  Qy  and  vice-versa. 

(ii)  Qy  is  a  smooth  semi-algebraic  submanifold  of  Q,  and  can  therefore  be 
stradfied  into  disdnct  fiber  bundles  [14].  By  a  bundle,  we  mean  the  set  of 
all  systems  which  differ  only  by  the  exact  values  of  their  eigenvalues;  but 
for  which  the  number  of  disdnct  eigenvalues  as  well  as  the  respecdve 


MODAL  STABILITY  OF  IMPERFECT  CYCLIC  SYSTEMS  U3 

orders  of  the  degenerate  eigenvalues  are  the  same.  Within  the  degenerate 
set.  Qy,,  the  crossing  from  one  bundle  to  another  can  also  lead  to  radical 
dynamical  changes.  Each  bundle  is  represented  by  a  specific  Jordan  block 
of  a  certain  order.  Note  that  each  bundle  is  also  transversal  to  Q. 


Theorem  2,2,  ref  [14], 

Qv  is  a  finite  union  of  smooth  sub-manifolds  with  codimension  satisfying 
Codim  Qv  >  2. 

Theorem  2.2  has  the  following  implications: 

(i)  Q/Qv  is  topologically  path  connected.  This  means  that  by  smooth  parame¬ 
ter  variations,  provided  that  the  number  of  variable  parameters  is  less  than 
the  codimension  of  2 v.  it  is  possible  to  smoothly  pass  from  one  member  of 
Q/Qv  to  another  without  reaching  any  singularity;  that  is,  without 
encountering  any  member  of  Qy.  Such  parameter  variations  will  typically 
not  lead  to  radical  dynamical  changes  in  Q/Qy. 

(ii)  Because  codim  Qy'>l.  it  follows  that  a  generic  one-parameter  family  of 
quadratic  systems  cannot  contain  any  degenerate  subfamilies.  Therefore 
under  one-parameter  deformadons  of  a  generic  family,  some  eigenvalue 
pairs  may  approach  each  other  but  cannot  be  coincident  (i.e.  cannot  col¬ 
lide).  After  approaching  each  other,  they  must  veer  away  rapidly,  giving 
rise  to  the  so-called  eigenvalue  loci- veering  phenomenon  [18],  under  one- 
parameter  deformadons  of  generic  families.  This  offers  a  theoredcal  expla- 
nadon  for  the  eigenloci  veering  phenomenon  which  has  been  observed  in 
perturbed  periodic  systems  without  a  corresponding  phenomenological 
base  [18.  19].  Furthermore,  this  phenomenon  holds  provided  the  system 
has  a  quadradc  structure,  irrespecdve  of  whether  the  model  arose  from  a 
condnuous  or  discrete  structural  system  [20]. 

This  rapid  eigenloci  veering  can,  under  the  right  condidons,  produce  the  mode 
localizadon  phenomenon  [18].  Since  the  dynamical  properties  of  any  linear 
constant-coefficient  system  are  totally  determined  by  its  eigen- structure  (eigen¬ 
values  and  eigenvectors),  and  since  the  eigenvalues  are  condnuous  funedons  of 
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Che  outrix  elements,  it  follows  that  radical  changes  in  the  dynamical  properties  of 
a  given  system  under  difierential  parameter  perturbations  ensue  principally  from 
a  large  disorientation  between  the  eigenvectors  of  the  tuned  (unperturbed)  and 
mistuned  (perturbed)  systems.  We  study,  in  Secdon  IV,  the  variadon  of  eigenvec¬ 
tors  of  generic  families  under  di&rendal  random  parameter  -^erturbadons. 


III.  Bounds  on  Amplitude  Ratios 

Consider,  again,  the  equadon  set  for  the  dynamics  of  quadradc  systems: 

Mxo*Kxo=f,  (3.1) 

where  M  and  K  are  posidve  definite  matrices.  For  tuned  cyclic  systems,  M  and  K 
have  the  addidonal  propeny  of  being  circulant.  Taking  the  Laplace  transform  of 
(3.2)  under  zero  inidal  condidons,  gives: 

^K)Xo{s)  =  F(s),  (3.2) 

or 

A(s)Xo(s)^Fis)  (3.3) 

where  4-  K.  Suppressing  s  in  all  subsequent  calculadons  leads  to: 

(3.4) 

The  posidve  definite  nature  of  M  and K  guarantees  that  A~^  exists  for  all  5  on  the 
Nyquist  contour.  Under  normal  operadons  of  the  system,  suppose  A  varies  to 
A  +  AA  :=  A,.  Let  Xq  then  change  to  )fo  •=  ^e-  Then,  for  the  same  excita- 
don  force  as  in  the  tuned  state. 

Xe  =  (A  +  AA)"‘  f.  (3.5) 

The  physical  nature  of  the  system  guarantees  chat  A  -t-  AA  will  always  remain 
symmetric  but  not  necessarily  circulant  since  a  true  mistuning  destroys  cyclicity. 
Equadon  (3.5)  can  be  rewritten  as: 

X,  =  (A  +  AA)"‘  •  f  =  (/  +  A~^AA)"‘  ■  A'^F.  (3.6) 


Subsdtudng  (3.4)  into  (3.6)  gives: 


modal  stability  of  imperfect  cyclic  systems 

X.  =  (/  -Xq. 

NormaUy  AAwUIbea  dififerentiai  perturbation  of  A.  so  that: 

P(A-'aA)<  1. 

where  p(-)  is  the  spectral  radius  of  (•).  Hence 

(/+4-‘^r‘=  £(-i)*(a-‘m)* 
Substituting  (3.8)  into  (3.7)  gives; 

•^«  =  f  ATo. 

Taking  norms  in  (3.9)  gives: 

IIX^W  =  II  £  (-1)'‘(A"‘aA)*A'oII 
i=0 
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(3.7) 


(3.8) 


(3.9) 


5  2  II  A'^AAII*-  llXnll 

i=0 


(3.10) 


Let  ))A"^;i4  II  =  r.  Because  AA  is  a  difierential 
thatr  <  1.  Hence: 


pemubation  of  A,  it  follows 


X*ltSII3(oll£  r*=  llj^oll  l+r  +  r2+  ...  ^ 
k^O 


IIXqII 

^  .  since  r  <  1. 


(3.11) 


1 ”  I- iia-'aa 


(3.12) 


Write 
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A  =zD  ■>■€  =DU  *  D-^C) 


(3.13) 


where  4.  is  a  diagonal  matrix  of  the  uncoupled  system  dynamic  matrix  and  C  is 
the  relative  coupling  dynamic  matrix,  such  that  the  minimum  eigenvalue  of 
D~^  C  at  any  ^equency  gives  the  coupling  index  of  the  system  at  that  frequency 
[21].  If  the  norms  in  (3.12)  are  //“-  norms,  then,  over  the  frequency  interval  Q: 


ess.sup 

U)£0 


;iX,(a)jii, 


^ess.sup 

oxa 


1 

<3„„AA(oi) 


(3.14) 


where  Orntx!  )  Omuil  )  correspond  to  maximum  and  minimum  singular  values 
of  (■)  respecdvely.  .Note  that  all  the  m-'mces  and  vectors  considered  above  are 
functions  of  frequency  r=i  a). 

Because  A  is  symmetric  it  follows  from  (3.13)  that: 

^min  (A)  =  o^(i)[/  +  D-'C]) 

=  KmiD)-KmiI  +D-^C), 

=  (3.15) 

by  the  eigenvalue  shift  theorem,  where  dmm  is  the  minimum  eigenvalue  of  D.  Let 

,  ilX,(co)ll, 

dnun  =  a,  and  =  Xq.  At  any  frequer'-v  co,  let  =  n(o)). 

llAo(<^)  ll- 

Then  (3.14)  reduces  to: 


5,  =  ess.sup  n(0))  2 

(oeQ 


1  -  ess.sup 
coeQ 


_I _ 

am«J^(a)) 
a((i))(l  +  /io(o))) 


(3.16) 


where  £2  is  a  frequency  interval  of  interest.  In  some  cases  it  is  possible  to  define 
£2  by  the  semi-open  interval  £2  =  [  0  "»).  Here  Xo  is  called  the  coupling  index  of 
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the  system.  The  system  is  decoupled  when  Xo  =  0.  It  is  weakly  coupled  if 
Xo  <  1.  and  is  strongly  coupled  if  Xo  2  1.  In  general,  0  2  Xo  ^  Observe  that 
Xo(k,  (0)  is  a  function  of  both  the  structural  coupling  k,  and  frequency  to.  Inequal¬ 
ity  (3.16)  leads  to  the  following  conclusions: 

(i)  The  misnined  to  tuned  amplitude  ratio  is  determined  by  the  maximum  peak 
of  the  mistuning  strength  an,„<14(a}).  the  minimum  strength  of  the  weakest 
link  in  the  system  aid),  and  the  minimum  peak  of  the  coupling  index 
(strength)  Xo(a)). 

(ii)  A  variation  in  rigidity  (coupling)  affects  the  ratio  of  (3.16)  monotonically 
for  fixed  values  of  a^„(,AA)  and  a.  This  is  because  at  any  given  fre¬ 
quency,  Xo  vanes  continuously  and  monotonically  as  the  coupling  is  varied 
[13]. 

(iii)  A  reduction  in  a  caused  by  a  reducdon  of  mass  of  the  blades,  and/or  more 
flexible  blades,  increases  the  ratio  (3.16)  monotonically.  More  specifically, 
at  any  frequency  when  Xq  ^  0,  from  (3.16): 


S(  ^  ess.sup 


_ _ 

a(a))-Om„AA(a)) 


>  1,  forOm„AA((o)  >0,  V  toefl. 


Hence  under  weak  coupling  across  the  frequency  interval,  the  amplitude 
ratio  depends  entirely  on  the  relationship  between  the  frequency  response 
of  the  mistuning  strength  and  that  of  the  strength  of  the  weakest  blade  in 
the  assembly.  Under  these  conditions,  the  maximum  amplitude  ratio  will 
arise  from  the  blade  with  the  worst  mistune  [22]. 

IV.  Eigenvector  Rotations 

In  section  II,  we  showed  that  generic  systems  QiQy  will  typically  have  dis¬ 
tinct  eigenvalues,  while  degenerate  systems  Qy  will  typically  have  repeated 
eigenvalues.  To  study  eigenvector  perturbations  for  generic  systems,  regular 
analytical  methods  will  work,  while  for  eigenvector  variations  in  the  system  Qy 
we  require  singular  perturbations  [23].  Let  A  e  C"*"  be  the  dynamic  matrix 
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arising  from  any  system  Q,s  QlQ^-  Let  F  represent  the  parameter  space  and  let 
|i  €  r  be  a  p-dimensionai  parameter  vector.  If  Codim  Cv  ^  ^  then  for  any 
^c  Fc  E''.  where  p<r,  differential  parameter  variations  in  A(^)  will  not  lead  to 
eigenvalue  degeneracies.  Thus,  if  the  eigenvalues  of  A(ix),  given  by 
Xi(n),  X2(ii),  •  •  •  ,X„(p.),  are  distinct  when  p.=0  they  will  continue  to  remain  dis- 
dnct  when  |i  is  small,  by  continuity  arguments. 

Let 


/l(6ii)  =  A(0)  +  &4, 


(4.1) 


where; 

5A  =5^*- A  (11)1^1.0  (4-2) 

5A  can  be  expanded  in  Taylor  series  form  as:  SA  =  AA  +  higher  terms  dependent 
on  p..  To  a  first  order  approximation  we  can  write  the  perturbed  matrix  as: 

A=Ao+AA.  (4.3) 


Write 


Ao  =  UA 

where  U  is  the  modal  matrix  of  Aq,  and  V'  =  U~^  where: 


(4.4) 


(/  =  («,.  ui, 


and 


V  =  [vi  V2  .  ■ • • 


with 


A  =  diag(Xi,  ^2,  •  •  • ,  A.,). 

where  (■)*  is  the  complex  conjugate  transpose  of  (■). 

Since  Aq  is  also  generic,  we  can  write  the  perturbed  modal  expression  as: 

Ao  +  AA  =  (C/  +  Af/][A  +  AA][C/  +  AC/r^  (4.5) 

where  AU  is  the  perturbation  in  U  resulting  from  AA  while  AA  is  the 
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conesponding  perturbation  in  A  resulting  from  AA.  Under  eigenvector  normali¬ 
zation,  II  u,  II  =  I,  and  II  u,  t  Aui  II  =  1.  Equation  (4.5)  can  be  solved  as: 

AqU  +  AAU  +  A  qAU  —  U\  +  U  AA  +  A£/  A  (4.6) 

where  we  neglect  second  order  terms  like  A£/AA  and  AAAU  [24],  As  a  measure 
of  the  eigenvector  vanadons,  we  begin  by  writing  Ak,  as  a  linear  combination  of 
all  the  eigenvectors  since  the  eigenvectors  span  the  whole  n-dimensional  space. 
Thus: 

lu,  =  £  Iji  Uj  (4.7) 

y=i 

Or: 


AU  =  UL.  (4.8) 

Now  solving  for  AA  in  (4.6)  gives: 

AA  =  AqU  +  AAU  Aq  AU  -  A  -  £/"'  Af/A.  (4.9) 
Observe  that  AqU  -  A  =  0,  so  that 

AA  =  AAU  +  AL-LA.  (4.10) 

Notice  that  the  diagonal  elements  of  (AL  -  LA)  are  zero.  Hence: 

AX,  =  [(/"^AAC/lii  =  v*  AAUj. 

To  solve  for  AU,  we  need  L.  The  off-diagonal  elements  of  L  are  given  by  Skelton 
[24] 

Iji  =  (kj  -X,)”‘  V,*  AAUj  for  i*j,  ij  =1,2  ■  ■  ■  n  , 


or: 

iij  =  (^1  “  Xy)"'  V*  AAUi  for  i*j,  i,j  =  1,2  ■  ■  ■  n. 

To  determine  observe  that  the  constraint  equation  II  Ui  +  Auj  11  =  1  contains 
Thus 


Hu,  +  Au,  II  =  (<  Uj  +  Aui,  u,  +  AUi  >  =  1. 


(4.12) 
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Or 


u*  Ui  +  2m*  Am,  Am,*  Am,  =  1 . 

But  Ui  Ui  ^  1,  so  chat 

2m,*  SUi  Am*  Am,  =  0. 


Therefore: 


n 

Am,  =  ^  ~  X  h* 

y  =  l  ;.l 

I 


=  X,  +  /„  Mi 


(4.13) 


(4.14) 


where: 


Thus: 


Letting: 


=  I 

y-i 
y  »* 


h‘  ‘*J- 


l}i  +  (2  +  2u*Xi)lu  +  (2m,* Xi  +  x,*Xi)  =  0 


uUi 


(4.15) 


(4.16) 


and 

yi  =xU,. 

gives  (on  accepting  the  positive  solution  of  the  quadratic): 

/u=-(l  +  zi)  +  (l  +  z?-yi)^  (4.17) 

Since  the  eigenvectors  Mi  and  u,-  -t-  AMi  can  be  nomalized  to  unity  and  since  each 
vector  is  represented  by  a  magnitude  nii  and  an  angle  6,,  the  natural  measure  of 
modal  variations  is  9,  since  miSl  after  normalization.  Knowing  ail  the  elements 
of  L,  we  can  now  determine  6i  as: 
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<Ui,  U,  +  Ai4,>  =  II  Ui  II  ll«,-  +  Alii  II  COS  0;.  (4.18) 

But  Uu,-  II  =  II  Ui  *  Aui  II  =  1. 

Hence: 


cos  0i  =  <Ui,  Ui>  +  <Ui,  Aui> 


=  [  +  u-  (Xi  +  /„  Ui) 


=  I  -r  u'x.  +  /,i 

=  (1 OS0iSJl/2.  (4.19) 


Consequently  for  the  occurrence  of  no  vector  rotation  under  parameter  variations, 
we  require: 

z?-yi=0  (4.20) 


Or: 


X-  Ui  u-Xi  - X-Xi  =  x' (UiU*  - /)Xi  =  0  (4.21) 

This  implies  Xi  belongs  to  the  null  space  of  (UiU*  -  /),  that  is: 

iUiU^ -I)i  IJiUj^O  (4.22) 

y  »  I 

where: 

Iji  =  (ki-\j)~^  V*  ^  Ui  i*j 

The  nearer  the  expression  (4.20)  is  to  zero,  the  less  the  corresponding  eigenvector 
rotadon  under  the  given  perturbadon.  Let 


oq  —  2  Iji  Uj  (UiUi  ~f)  2  ^y*  i  ~  1,  2  ■  ■  ■  n. 

y  ■  <  y  - ' 

i*i  j0t 

Then  max(ai )  gives  the  eigenvector  with  maximum  rotadon. 

i 

The  conclusions  are  the  following: 
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a)  If  the  separation  between  the  eigenvalues  is  very  large,  (i.e.  (X;  -  Xy)  is 
very  large  for  all  i,j),  then  Iji  is  relatively  small  ai,J  eigenvector  roution 
will  be  correspondingly  small. 

b)  If  v*  AAu,'  =  0.  then  eigenvector  roution  will  also  be  relatively  small,  pro¬ 
vided  Iji  ^  <">. 

For  example,  if  Aq  is  Hermidan  as  is  the  case  in  all  quadratic  systems,  and 
M  =  a/,  ae  C,  then 

V*  ^Au,  s  0,  V  ij  =  1,2,  ■  ■  •  n. 

Thus,  identical  increases  or  decreases  in  the  diagonal  elements  of  a  quadradc  sys¬ 
tem  will  not  produce  unexpected  amplitude  excursions  [25]  because  it  caiuot 
produce  eigenvalue  splitdngs  in  formerly  degenerate  families.  Therefore,  such 
perturbadon  cannot  take  a  system  either  across  the  boundary  of  the  bifurcadon  set 
or  across  diSerent  bundles  of  Q^.  Geometrically,  this  implies  that  degenerate 
eigenvalues  in  systems  belonging  to  a  bundle  in  Qv  cannot  be  lifted  by  perturba- 
dons  that  leave  the  perturbed  system  in  the  same  bundle  of  Qy.  Indeed,  define 
the  eigenvector  sensidvity  matrix  of  a  quadradc  system  as 

S  =  =  L, 

from  eqn.  (4.8).  Defining  the  eigenvector  sensidvity  metric  measure  by 

=  ll5ili=  i  Iji 

<  >  i 

where  Sp  is  the  Frobenius  norm  of  5  shows  that  the  maximum  eigenvector  sensi¬ 
dvity  is  obtained  at  the  posidons  of  minimum  eigenvalue  separadon,  which  is  not 
difficult  to  compute.  Altemadvely,  (5/r)max  ^so  corresponds  to  the  posidon  of 
maximum  angular  rotadon  between  the  tuned  and  mistuned  system  eigenvectors. 
This  condidon  is  evidenced  by  strong  eigenloci  deformadons. 

If  A  (CD)  is  a  frequency  response  matrix  arising  from  a  generic  system,  the 
eigenvalues  Xi(o))  and  eigenvectors  UiCco)  are  also  condnuous  funedons  of  fre¬ 
quency.  We  can  therefore  plot  the  frequency  response  funedons  5/r(co)  to  deter¬ 
mine  the  fir^'Mencies  at  which  maximum  deformadons  take  place. 
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V.  ExamplM 

To  illustrate  the  theory  so  far  developed  we  consider  two  examples.  The  first 
is  an  interconnected  linear  chain  of  oscillators.  This  has  been  studied  by  Arnold 
[13]  and  more  recently  by  Piene  [18). 

Example  1:  fvtode  Localization  in  Generic  Periodic  Systems. 

Consider  a  coupled  pendulum,  as  shown  in  Fig  I,  with  identical  masses  but  of 
different  lengths  /i  and  li,  where  li  is  a  perturbation  of  /j,  i.e.,  /2  =  +  ^i)-  If 

we  put  =  /,  then  the  kinedc  energy  is  given  by 

r  =  7  9i]  (5.1) 

while  the  potential  energy  is  given  by 

0  =  m/~  +  m(/ A/)~— +  7(0j  -  02)^.  (5.2) 

«•  A  ^ 


Fig.  1.  Two  coupled  oscillators. 


Under  unit  gravitational  force,  application  of  Lagrange’s  equations  results  in  the 
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equation  of  motion; 


where 


Vf 


and 


^  = 


■fll 

■01 ' 

02 

-  [K] 

02 

ml^ 

0 

0  mil 

+ 

ml  *  k 

-k 

=  0 


-k  m(l  +  k 


The  dynamic  matrix  for  tne  above  system  is  given  by 
A(0})= 


ml  *  k  -  ~k 

-k  mil  +  d/)  +  k-mXil  -h 


Rewrite  A  (at)  as: 


Ai<a) 


=  [-“ 


a  -k 
k  b 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


(5.7) 


The  characteristic  equation  of  A  (co)  is  given  by 

X^-(a*b)X  +  iab~kh::0 

Both  Af  and  K  are  symmetric  and  positive  definite.  The  eigenvalues  of  A  (  )  art: 


,  f(a  +  b)±  -  b)^  1 

^1.2  “ 


(5.8) 


Note  that  Xi,2  cannot  be  degenerate.  Thus  under  one-parameter  deformations,  the 
eigenvalues  can  deform  but  cannot  collide. 

Indeed, 


do  2 


I 


a-b 


V{fl  -b)^  +  Ak^  J 


(5.9) 


and 
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Hence: 


V(a  -6)2  + 

da  da 


da  2 


da 


1 

2 


when  a=b. 


(5.10) 


(5.11) 


The  distance  between  the  eigenvalues  is  given  by: 

‘^x=!Xi-X2l=  V(a-6)2  +  4ik2  (5.12) 

which  assumes  its  minimum  value  of  2k  when  a-b  or  when  |dX,/da|  =  Vi  This 

represents  the  tuned  state  of  the  linear  chain.  For  a  fixed  mistuning  value  {a-b), 
d\  depends  essennally  on  k.  If  {a-b)  is  small,  it  is  clear  that  5f  -♦  <*  as  X  -♦  0. 
The  modal  matrix  of  the  chain  is  given  by 


-1  -1 

(g  -  6)  -  V(a  -  b)^  4)^  (a-b)  +  V(a  -b)^  +  4k^ 
2k  Ik 


(5.13) 


O:  serve  that  u]  UjsO,  V  k,a,b.  Under  tuned  conditions,  a=b,  then 

V.  =  [:1 1']  (5.14) 


However,  consider  the  very  interesting  situation  when  the  mistuning  to  coupling 
ratio  is  rather  large.  That  is  to  say; 

^^’1. 

k 

Then  (a-b)^  »k^,  and  k^  becomes  negligible  in  the  eigenvector  expressions. 
Expanding  the  term  under  the  radicals  and  neglecting  second  and  higher  order 
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terms  gives: 

12 

=  (a-6)[l  >  7  — ^  +  •  ••  (5.15) 

2  a  -  b 


In  this  case  the  modal  matnx  reduces  to: 

-1  -1 

k  a  -  b 
a  -  b  k 

An  energy  exchange  now  cakes  place.  The  second  component  of  the  first  mode 
becomes  vanishingly  small  while  the  corresponding  component  of  the  second 
mode  becomes  extremely  large.  This  is  an  extreme  case  of  classical  vibration 
absorber,  and  is  the  mode  localization  phenomenon.  We  therefore  conclude  that 
mode  localizadon  (or  extreme  energy  exchange)  will  occur  in  a  generic  system 
under  one-parameter  deformations  if  the  following  conditions  are  satisfied: 

•  at  any  frequency  cu  where  the  system  is  almost  decoupled,  i.e.,  X(a))o  -*  0. 
(Note  that  X(0))o  -♦  0  as  k-»0). 

•  when  the  mistuning  to  coupling  ratio  ^*7^  >  1. 


(5.16) 


At  the  localization  stage  the  eigenvalue  and  eigenvector  sensitivities  take  on  their 
maximum  values,  i.e.  both  IIAA  A~^  11^  and  II  Af/  U~^  11^  have  their  maximum 
values.  Localized  modes  always  produce: 


(5.17) 


Example  2:  Cyclic  Systems. 

Grnsider  three  identical  masses,  m,  arranged  in  a  ring  structure  and  intercon¬ 
nected  by  identical  springs  kg.  Assume  that  ail  the  masses  are  hinged  to  the 
ground  by  torsional  springs  of  strength  k„  and  that  the  radius  of  the  ring  is  r ;  as 
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shown  in  Fig.  2. 

The  basic  equations  of  motion  of  this  '‘ring”  is 

S4x  *  Kx  =  f 


where 


2kt  +  — 

-kc 

-kc 

m 

0 

0 

kt 

0 

m 

0 

.  K  = 

-*c 

Ik,^  — 

-kc 

0 

0 

m 

r 

-kc 

-  kc 

kf 

(5.18) 


(5.19) 


Fig.  2.  Models  of  (a)  the  cyclic  chain,  (b)  the  linear  chain 
with  three  degress  of  freedom. 
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Using  group  theoretic  arguments  [13],  we  can  easily  deduce  that  the  above  $ys> 
tern  has  degenerate  eigenvalues  occurring  as  doublets,  by  cyclicity  of  the 
corresponding  system  matrices.  Consequently,  every  quadratic  cyclic  system 
Qe  c  fiv  Furthermore  all  perturbations  of  the  above  system  preserving  the 
cyclic  structure,  leaves  the  modal  geometry  invariant  [3,  25].  Indeed  the  eigen¬ 
values  of  the  above  system  are  given  as: 


\  —  A-  —  +  —  X-  —  +  — 

*  w  ’  w  m  '  ^  mr  m 


(5.20) 


Write  the  dynamic  matrix  of  this  system  as: 

where 


a  -b  -b 
-b  a  -b 
-b  -b  a 


a  =  2kc  k,  -  oi"m,  b  =  kc- 
We  now  consider  a  diagonal  perturbation  of  the  form: 

£  =  diag  (ei,  ei,  ^3). 

Then 

=  Ao((0)  +  E. 


(5.21) 


(5.22) 


(5.23) 


This  would  correspond  to  the  realistic  situation  where  there  are  slight  changes  in 
the  values  of  the  ground  spring  I;,,  depending  for  example  on  how  the  blades  are 
coupled  to  the  disk  in  bladed  disk  assemblies  [22].  The  major  difference  between 
the  behavior  of  (degenerate)  cyclic  systems  and  generic  systems  are  the  follow¬ 
ing: 

(i)  For  generic  systems,  all  the  eigenvalues  and  the  distance  between  adjacent 
pairs  inaeases  as  the  coupling  kg  increases.  Consequently  the  probability 
of  mode  localization  decreases  as  kg  and  hence  >.((0)0  increases.  On  the 
other  hand,  perturbations  which  split  the  degenerate  eigenvalue  of  cyclic 
systems  turn  them  into  generic  systems  with  pathologically  close  eigen¬ 
values  [7].  Hence  for  previously  cyclic  systems  whose  eigenvalues  bifur¬ 
cate  under  perturbations,  Sf  is  very  large.  Therefore  such  systems  are 
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susceptible  to  mode  localization,  independent  of  the  values  of  the  coupling 
strength /lie.  Recall  that 

S/r  =  IILII^. 

where 

li!  =  0^!  -  AAu^.  i^j.  i,j  =  1,  2,  ■  ■  ■  n. 

(ii)  Consequently  the  only  way  to  avoid  large  values  of  S/r  in  such  a  situation  is 
if  and  only  if  II  v*  AAUj  II  as  0  or  in  the  neighborhood  of  zero.  Perturbations 
that  induce  this  condidon  are  precisely  those  that  will  not  induce  radical 
dynamical  changes  in  misnined  cyclic  systems.  It  was  already  shown  that 
if  dA  =  a/  ,  then  II  v'  SAu^  II  s  0 

(iii)  Of  the  remaining  possible  perturbadons  those  that  have 

II V*  AAUj  II  =  e  <<  I  will  produce  minimum  dynamical  changes.  All  oth¬ 
ers  for  which  II  v  *  II  is  not  small  will  give  suscepdbility  to  mode  local- 

izadon,  no  matter  how  strong  the  interblade  coupling. 

The  following  numerical  example  amplifies  the  above  observadons.  We  con¬ 
sider  the  case  of  the  so-called  ‘strong  coupling’,  using  the  following  values: 
kf.=9.5,  k,-l,  a  =  20,b~  9.5,  63=0,  Clearly,  the  rado  of  mis- 

tuning  to  coupling  strength  is  very  small.  Now,  in  order  to  compute  the  frequency 
response  curves,  we  need  some  damping  to  obtain  finite  amplitudes  at  resonance. 
Assume  hystetedc  damping  of  0.01  for  all  cases.  Without  loss  of  generality,  the 
response  to  be  computed  is  the  direct  receptance,  i.e.  the  response  of  each  node  to 
individual  excitadon.  We  turn  the  ring  into  a  linear  chain  by  putting 
b  =  ki2  =  ^31  =  0  in  equadon  (3.21).  Then  Aq  becomes  a  tridiagonal  banded 
matrix. 

The  frequency  response  of  the  tuned  and  mistuned  systems  of  the  linear  chain 
are  shown  in  Fig  3.  The  illustradon  is  windowed  around  one  of  the  resonant  fre¬ 
quencies  of  the  coupled  system.  Nonce  that,  at  the  tuned  state,  ..le  amplitudes  of 
nodes  1  and  3  are  equal  on  account  of  symmetry,  while  that  of  node  2  is  double 
that  magnitude. 

Because  the  system  is  now  generic,  and  therefore  exhibits  modal  stability,  all 
nodes  have  almost  the  same  response  panems  and  magnitudes  as  in  the  tuned 
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system.  This  is  also  the  case  when  we  change  the  sign  of  e  2 .  from  -0. 1  to  0. 1 . 

When  we  repeat  exactly  the  same  procedure  for  the  circulant  system,  a  very 
different  picture  is  obtained.  Fig  4  shows  the  response  of  individual  nodes  com¬ 
pared  with  the  tuned  case.  This  case  corresponds  to  a  2-parameter  perturbation, 
with  Cl  =0.1,e2  =  ~01.«3  =0. 


(d)  nodtX 


(b)  node!  (c)  nodel 


Fig,  3,  Effect  of  mistuning  on  ih^  response  curves  of  the  linear  chain.  Note  the 
preservation  of  the  shape  of  the  curves  around  resonance,  and  the  minimal 
difference  in  the  peak  amplitudes  of  the  tuned  and  mistuned  systems  (■•  -  -  tuned 
systems; _ mistuned  system). 
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(a )  node  I 


Fig.  4.  Effect  of  two  parameter  mistuning  on  the  response  curve  of  the  cyclic 
chain.  Note  the  severe  reduction  in  the  amplitude  at  node  3,  which  is  only  50%  of 
the  tuned  system  ( - tuned  systems; _ mistuned  system). 


Notice  that  the  node  with  zero  toistunittg  (mode  3)  now  has  a  reduction  in 
amplitude  of  almost  50%.  This  extremely  unequal  amplitude  distortion  (Fig  4)  is 
the  case  no  matter  how  small  the  magnitude  of  the  perturbation  is,  so  long  as  we 
keep  the /arm  of  mistuning,  and  the  mistuning  does  not  actually  vanish. 

If  we  now  change  the  mistuning  matrix  in  a  very  small  way,  by  making 
€2^0.1,  we  obtain  the  response  curves  in  Fig  S.  We  now  notice  a  substantial 
difference  in  the  geometry  of  the  curves  in  Fig  5,  compared  to  those  in  Fig  4. 
Thus,  a  very  small  change  in  the  perturbation  matrix,  now  results  in  a  consider¬ 
able  difference  in  the  vibration  response  at  the  individual  nodes.  The  question  of 
which  node  will  be  most  responding,  or  the  one  having  the  least  amplinide,  is  now 
not  as  easy  as  one  would  have  expected.  In  Fig  4,  it  is  node  3,  while  it  is  node  1  in 
Fig  S.  In  fact,  the  amplinide  of  node  3  has  been  increased  by  about  1(X)%  from 
Fig  4  to  Fig  5,  merely  by  changing  only  one  entry  in  the  system  matrix  from  19.9 
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to  20.1,  a  change  of  less  than  1%  ! 

The  foregoing  examples,  based  on  a  simple  3  degrees  of  freedom  model  of  a 
circular  ring  or  disk  only,  illustrates  the  instability  induced  by  cyclicity.  It  is 
clear  that  the  qualitative  conclusions  to  be  drawn  from  Fig  4  are  inconsistent  with 
those  from  Fig  5,  although  the  difference  between  the  two  mistuned  matrices  is 
very  small  indeed.  We  emphasize  that  these  results,  obtained  for  just  a  cyclic 
chain,  are  not  necessanly  applicable  to  bladed  disks  in  all  generality,  especially 
those  models  in  which  cyclicity  is  ignored.  However,  when  bladed  disk  systems 
are  well- modeled  to  include  the  effects  of  blade  coupling,  blade  or  isk  mistuning 
and  cyclicity,  similar  distortions  in  the  geometry  of  the  frequency  response  curves 
can  result.  The  subject  is  currently  under  investigation  by  us. 


(d)  nod€  I 


Fig.  5.  Effect  of  one-parameter  mistiming  on  the  response  curve  of  the  cyclic 

chcan.  Note  the  symmetrical  unfolding  of  the  degenerate  singularity  ( - tuned 

systems; _ mistuned  system). 
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VI.  ConeitMions 

(i)  For  generic  systems,  to  which  linear  periodic  chains  of  oscillators  belong, 
differential  parameter  perturbations  are  significant  for  the  system  dynamics 
only  under  weak  coupling  conditions  when  the  mistuning  to  coupling  ratio 
exceeds  unity  (Example  1).  Under  all  other  conditions  that  do  not  induce 
eigenvalue  degeneracy;  small  magnitudes  of  mistuning,  or  the  type  of  mis¬ 
tuning,  is  irrelevant  to  system  dynamics. 

(ii)  For  degenerate  systems  to  which  a  tuned  cyclic  system  with  circulant 
dynamic  matrices  belongs,  it  is  not  just  the  mistuning  to  coupling  ratio 
which  is  significant  in  the  determination  of  the  perturbed  system  dynamics. 
The  type  of  mistuning  assumes  a  far  greater  importance  than  the  mistuning 
to  coupling  ratio.  All  types  of  mistuning  that  move  the  system  either  across 
the  boundary  of  the  bifurcation  set,  or  from  one  fiber  bundle  of  the  degen¬ 
erate  set  to  another  within  Qy,  will  lead  to  topological  catastrophes  [15]. 
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APPENDIX  4 


The  Stability  of  Frequency  Response  Curves 


SUMMARY 


In  this  Appendix,  we  highlight  one  of  the  results  obtained  so  far,  namely:  that  the 
modes  of  vibration  of  qrdic  structures  are  imstable  under  arbitrarily  small 
perturbation.  This  instability  is  not  the  usual  ill-conditioned  problem  of 
numerical  analysis.  It  has  nothing  to  do  with  the  computational  algorithm.  The 
eigenvector  instability  results  because  the  perfect  cyclic  system  is  "physically  ill- 
conditioned",  since  a  very  small  perturbation  changes  its  dynamics 
characteristics  dramatically.  This  is  significant  because  many  aerospace 
structures  have  circular  profiles.  The  implication  of  eigenvector  instability  for 
modal  control,  forced  response  amplitudes,  sensitivity  analysis,  etc,  therefore 
needs  further  investigation. 


Eigenvector  Stability,  Forced  Response,  and  Turbine  Blade  Failure 

The  structural  integrity  of  turbine  blades  used  in  jet  propulsion  systems  is 
sometimes  compromised  by  the  rare,  but  very  dangerous,  failure  of  some  "rogue 
blades".  This  problem  has  been  addressed  by  different  investigators  of  the 
mistuning  problem.  However,  they  often  obtained  conflicting  results. 

This  is  because  the  response  obtained  from  each  rotor  studied  by  each  author 
depenck  on  the  eigenvectors  of  the  rotor  system  matrix  A.  In  general,  the  matrix 
A  will  be  different  for  each  model  used  by  each  author,  although  the  difference 
may  be  very  small.  In  fact,  mistuning  is  usually  small. 

However,  the  problem  created  by  mistuning  is  not  always  small.  Thus,  although 
each  A  in  the  family  is  difierentiably  dependent  on  the  mistuning  parameter  e  in 
the  neighborhood  of  the  origin  of  E,  the  corresponding  eigenvectors  is  not. 
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Consequently  two  almost  similar  rotors  may  produce  dramatically  different 
vibration  responses,  if  their  respective  system  matrices  are  different 
perturbations  of  the  same  nominal  matrix.  An  effective  demonstration  of 
unstable  frequency  response  curves  in  a  simple  3  degree  of  freedom  cyclic 
system  is  given  in  the  following  examples.  First,  we  examine  the  instability  of 
eigenvectors,  then  the  instability  of  frequency  response  curves. 


Numerical  Example;:  illustrating  Eigenvector  Instability 


Example  1 

At  least  three  coordinates  are  required  to  define  a  cyclic  system  uniquely. 
Therefore,  we  consider  the  simplest  possible  example:  a  3x3  drculant  matrix 
with  real  elements,  a,  be  R. 


aO  = 


a  b  b' 
b  a  b 
b  b  a 


(A4.1) 


Using  the  following  perturbation  matrices,  where  se  R  is  a  very  small  parame¬ 
ter,  we  can  generate  two  matrices  A\  =<i0  +  £i  and  A2*fl0  +  £2  that  are  very 
close,  and  such  that  these  depend  smoothly  on  e,  and  ase-^O,  ->flO  <-  Ai. 
Thus,  if 


then 


'  e 

0  O' 

e  0  0 

£i  = 

0 

0 

e  0 
0  0. 

£2  = 

0  -£  0 
0  0  0 

(A4.2,A4.2) 


a  b  b 
b  a  +e  b 
b  b  a 


(A4.4) 


and 


^2  = 


a  +e  b  b 
b  a-e  b 
b  b  a 


(A4.5) 
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Note  that  I  I  i4i  I  I  =  11^211,  where  1 1  •  I  I  is  some  norm. 

Consider,  for  an  illustration,  a  situation  when  a  =  100,  b  =  45,  e  =  0.1.  We  can  com¬ 
pute  the  eigenvalues  of  A  i  and  A  2  respectively  as 

Ai  =  diag  (10.0667,  145.0333,  145.1000)  (A4.6) 

and 

A2  =  diag  (10.0000,  144.9423,  145.0578)  (A4.7) 

Notice  that  the  eigenvalues  of  the  two  matrices  are  very  close.  If  we  now  com¬ 
pute  the  corresponding  eigenvectors,  we  get 

■  .9993  -.5004  1.0000' 

1/1=  .9993  -.5004  -1.0000  (A4.8) 

[1.0000  1.0000  0 

■  .9985  -.2681  l.OOOO 

U2  =  1.0000  1.0000  -  .2678  (A4.9) 

.9993  -.7328  -.7313, 

We  now  notice  a  significant  difference  between  the  eigenvectors  at  modes  2  and  3 
of  matrices  A 1  and  A  2  respectively.  For  example,  there  is  no  node  (a  point  where 
displacement  is  zero)  in  the  third  mode  of  U 2,  whereas  there  exists  such  a  node 
in  Lfi- 


Example  2 

In  the  second  example,  we  consider  the  following  drculant;  its  elements  are  com¬ 
plex  but  the  matrix  is  not  symmetric.  It  may  be  regarded  as  a  deformation  of  a 
symmetric  drculant. 

■  200  +  i(-lO)  -  95  +  i(  -  5)  -  95  +  i(5)  ' 

A,=  -95  +  t(5)  200  +  i(-10)  -95  +  i(-5)  (A4.10) 

[  -  95  -t-  j(  -  5)  -  95  +  i(5)  200  +  i(  - 10). 

We  test  for  modal  stability  by  computing  the  eigenvalues 
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Ai  =  diag  10  + 1(  - 10),  2863398  +  i(  - 10),  303.6603  +  i(  - 10)  (A4.11) 


and  the  eigenvectors 


1  +  i(0)  -  .5  +  i(.866)  -  .5  +  i(  -  .866)' 

L/i  =  1  +  i(0)  1  + 1(0)  1  +  i(0)  (A4.12) 

1  +  1(0)  -  .5  +  i(  -  .866)  -  .5  +  i(.866) 


Now,  we  apply  a  very  small  perhirbation  to  the  matrix  A  i  to  get: 


‘201  +  i(  - 10)  -  95  +  i(  -  5)  -  95  +  i(5)  ‘ 

A2=  -95  +  1(5)  199  +  i(-10)  -95  +  i(-5)  (A4.13) 

-  95  +  i(  -  5)  -  95  +  i(5)  200  +  i(  - 10) 


It  is  clear  that  the  matrices  A  i  and  A  2  are  'close',  since  I  I  £  I  I  =0,  where 


-1+1(0)  0  +  i(0)  0  +  i(0)' 

E=Ai-A2=  0  +  i(0)  l  +  i(0)  0  +  i(0)  (A4.14) 

0  +  i(0)  0  +  i(0)  0  +  i(0) 


The  computed  eigenvalues  of  A  2  are: 


A2  =  diag  9.9977  +  i(  - 10),  286.3218  +  i{  - 10),  303.6807  +  i(  - 10)  kA4.15) 


Now,  notice  what  happens  to  the  third  eigenvector  of  A|  (eq.  (A4.12)),  as  a  very 
small  change  is  made  using  E,  eq.  (A4.13),  to  transform  it  to  A 2-  The  eigenvector 
matrix  of  A  2  is: 


■  0.993  +  i(O.OOO)  -  .474  +  i(0.820)  -  0.556  +  i(  -  .867)' 

1/2=  1.000  +  1(0.000)  1.000 +  i(- 0.00)  -  .444  +  i(-.861)  (A4.16) 

0.997 +i(-0.000)  -.531  +  1(0.817)  1.000  +  J (0.000) 


Again,  it  should  be  noted  that  a  very  small  change  in  the  matrix  Ai  induces  a 
sigiuficant  qualitative  difference  in  the  eigenvector  at  certain  modes  of  A^  (eq. 
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A4.10)  compared  with  the  corresponding  eigenvectors  of  A2,  (eq.  A4.16). 


From  Modal  Analysis,  sometimes  known  as  eigenfunction  expansion,  we  know 
that  the  forced  response  amplitudes  are  related  to  eigenvectors.  Thus,  if  the 
eigenvectors  are  unstable  under  arbitrary  perturbation,  then,  the  forced  response 
curves  will  also  be  unstable  imder  arbitrary  perturbation.  This  is  illustrated  in 
Figs  A4.1  to  A4.2  below. 


node  1 


node  2  node  3 


Fic  44. 1  Effect  of  onc-panmcicf  mistuning  on  the  response  curve  of  the  cyclic  mem. 
braii^Note  the  symmetrical  unfolding  of  the  degenerate  singularity. 

( - tuned  system:  _ misluned  system) 


FigA4>2  Effect  of  two-parameter  mistuning  on  the  response  curve  of  the  cyclk  nacm- 
'Bhuie.'  Note  the  severe  reduction  in  the  amplitude  at  node  3,  which  is  only  50%  of  the 


tuned  system 


( - tuned  system:  _ mistuned  system) 


